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ABSTRACT 


I 

Fully  developed  oil  flow  in  a  smooth  circular  pipe  at  a 
mean  Reynolds  number  of  about  2100  was  subjected  to  a  nominally  sinu¬ 
soidal  flow  modulation  at  frequencies  ranging  from  0.05  -  1.75  Hz. 

It  was  observed  that  flow  oscillation  increased  the  critical  Reynolds 
number  and,  under  certain  conditions,  even  brought  about  laminariza- 
tion  of  the  flow,  which  would  be  intermittently  turbulent  at  the  mean 
Reynolds  number  under  quasi-steady  (infinitely  small  oscillation  fre¬ 
quency)  conditions.  The  occurrence  and  extent  of  laminarization  was, 
however,  found  to  depend  on  factors  such  as  the  intermittency  of  turbulent 
puffs  in  the  mean  quasi-steady  flow,  frequency  of  oscillation,  etc.  Two 
series  of  experiments  were  performed.  In  one  series,  the  oscillatory  flow 
was  almost  completely  laminarized.  In  the  other  series,  the  oscillatory 
flow  was  fully  turbulent.  In  both  the  cases,  instantaneous  velocities  in 
the  flow  were  measured  using  Laser  Doppler  anemometry.  The  instantaneous 
velocity  was  decomposed  into  time-mean,  periodic  and  random  components 
employing  ensemble  averaging  techniques.  The  experiments  indicated  that 
the  laminarized  oscillatory  flow  behaves  very  similar  to  laminar  oscil¬ 
latory  flow  at  either  end  of  the  Strouhal  number  range  studied.  The  os¬ 
cillatory  turbulent  flow  was  found  to  depend  on  both  the  Strouhal  number 
and  the  ratio  of  the  oscillation  frequency  (f)  to  some  characteristic  fre¬ 
quency  Cft)  of  turbulence  in  the  flow.  The  design  of  the  present  experi¬ 
mental  facility  made  it  possible  to  stu<iy'--the  flow  at  f/f  *  1,  a  condition 
that  could  not  be  attained  in  most  previous  investigations.  It  was  found, 
that  at  this  frequency  of  oscillation,  the  Reynolds  stresses  generally  re¬ 
mained  frozen  at  an  average  state  during  the  entire  oscillation  cycle. 

The  turbulent  structure  showed  significant  departures  from  equilibrium  at 
all  times  during  the  oscillation  cycle.  As  a  result,  there  was  a  net 
change  in  the  time-mean  velocity  profile  near  the  wall  and  a  net  increase 
in  the  time-mean  wall  shear  stress  and  power  loss  due  to  friction.  The 
data  also  indicated  that  the  direct  interaction  between  oscillation  and 
the  turbulent  structure  was  essentially  confined  to  the  Stokes  layer. 


ix 


The  study  suggests  that  quasi-steady  turbulence  models  may  not  be  ade¬ 
quate  to  describe  unsteady  flows  when  the  time  scale  of  unsteadiness  is 
comparable  to  that  of  dominant  turbulent  eddies. 
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I.  INTRODUCTION 

The  study  of  unsteady  shear  flows  is  relevant  to  many  areas 
of  application  such  as  aerodynamics,  ship  hydrodynamics,  biofluid  mechanics 
and  wind  engineering.  Much  of  the  study  reported  in  the  literature  on 
unsteady  flows,  however,  concerns  laminar  flows.  Exact  solutions  are  avail¬ 
able  for  relatively  simple  unsteady  laminar  flow  situations  in  the  clas¬ 
sical  literature  [(see  Rosenhead  1963)].  Laminar  flows  in  pipes  due  to 
periodic  pressure  gradients  have  been  studied  by  Richardson  and  Tyler 
(1929),  and  Uchida  (1956)  and  others.  The  situation  with  regard  to  unsteady 
transitional  and  turbulent  flows  is,  however,  less  satisfactory.  While 
there  have  been  some  attempts  to  analyze  the  problem  of  stability  in  per¬ 
iodic  flows  (e.g.,  Kerczek  and  Davis  1974,  1975,  Hino  and  Sawamoto  1975, 

Hino  et  al  1976) ,  much  remains  to  be  studied  in  this  area.  Of  the  labor¬ 
atory  experiments  on  periodic  pipe  channel  flow  that  have  been  reported  in  the  lit¬ 
erature,  most  pertain  to  high  Reynolds  number  turbulent  flows.  The  ear¬ 
liest  of  such  experiments  were  reported  by  Schulz-Grunow  (1940).  More  re¬ 
cent  experiments  in  this  category  are  those  performed  by  Hirose  and  Oka 
(1969),  Lu  et  al,  (1973),  and  Acharya  and  Reynolds  (1975).  In  all  these 
experiments  fully  developed  turbulent  pipe  flow  was  perturbed  by  imposing 
a  periodic  modulation  in  discharge  at  a  prescribed  frequency.  Of  these, 
the  work  of  Acharya  and  Reynolds  involved  the  most  detailed  measurements 
such  as  Reynolds  shear  stress,  though  the  amplitude  of  flow  modulation  in 
their  experiments  was  less  than  5  per  cent  of  the  mean.  These  experi¬ 
ments  suggested  that  even  at  that  amplitude,  the  structure  of  turbulence 
can  be  significantly  affected  in  the  Stokes  layer  near  the  wall.  Their 
measurements  in  this  layer  were  very  limited  because  of  the  thinness  of 
the  layer.  Lu  et  al  experienced  considerable  difficulty  in  the  use  of  hot 
film  anemometry  as  well  as  in  analog  data  processing  in  their  experiments 
on  oscillatory  flow  of  water  in  a  circular  pipe.  Their  measurements  of 
mean  and  turbulent  velocities  did  not  lead  to  any  significant  conclusions. 

One  would  intuitively  expect  that  the  effects  of  unsteadiness  on 
the  flow  structure  will  be  stronger  on  flows  in  the  neighborhood  of  transi¬ 
tion  than  on  flows  at  very  large  Reynolds  numbers.  This  is  because 
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transition  process  can  be  sensitive  to  the  strong  acceleration/deceleration 
in  the  unsteady  flow.  The  experience  with  steady  flows  subjected  to  spa¬ 
tial  pressure  gradients  has  confirmed  that  pressure  gradients  have  a  signifi¬ 
cant  effect  on  the  critical  Reynolds  number.  Strong  effects  of  periodic 
flow  modulation  on  the  flow  characteristics  at  transitional  Reynolds  numbers 
have,  indeed,  been  observed  in  the  experiments  of  Sarpkaya  (1966)  and  in 
the  more  recent  studies  of  Gerrard  (197  L)  and  Hino  and  Sawamoto  (1975). 

From  a  study  of  the  growth  of  disturbances  in  sinusoidally  modulated  pipe 
flow  in  the  mean  flow  Reynolds  number  range  of  2000-5000,  Sarpkaya  conclu¬ 
ded  that  flow  pulsation  increases  the  critical  Reynolds  number.  Gerrard 's 
experiments  related  to  periodic  pipe  flow  at  a  mean  Reynolds  number  of  3700 
while  the  experiments  of  Hino  and  Sawamoto  pertaineed  to  purely  oscillatory 
pipe  flow  (i.e.,  at  a  mean  Reynolds  number  of  zero).  These  latter  two 
studies  were  mainly  qualitative  but  they  indicated  that,  in  general,  tur¬ 
bulence  is  inhibited  during  the  acceleration  part  of  the  modulation  cycle  and 
enhanced  during  the  retardation  part.  In  fact,  the  flow,  under  certain 
circumstances,  appeared  to  remain  laminar  during  the  acceleration  part  of 
the  cycle.  Indications  of  inhibition  of  turbulence  or  delayed  transition 
have  also  been  observed  in  some  of  the  in-vivo  studies  of  pulsatile  blood 
flow  in  the  aorta  of  mamals  (e.g.,  Falsetti  et  al  1972,  Kiser  et  al  1976). 
However,  the  instrumentation  and  data  reduction  techniques  vised  in  the  blood 
flow  studies,  as  well  as  in  the  other  laboratory  studies  mentioned  above, 
were  not  sufficiently  sophisticated  to  yield  detailed  quantitative  informa¬ 
tion  on  the  ..low  structure. 

There  is  one  class  of  unsteady  turbulent  flows,  however,  that 
has  been  studied  in  fair  detail  using  moderate  to  highly  sophisticated  data 
reduction  techniques.  This  is  the  unsteady  turbulent  boundary  layer  in  a 
periodic  free  stream.  Most  of  these  studies  (Karlsson  1959,  Houdeville  et 
al  1976,  Cousteix  et  al  1977,  and  Patel  1977)  pertain  to  zero  pressure  gra¬ 
dient  boundary  layers  though  some  experiments  on  adverse  pressure  gradient 
flows  have  been  reported  recently  (Kenis>n  1977.  Houdeville  and  Cousteix 
1978).  The  experiments  on  the  zero  pressure  gradient  boundary  layers  have, 
generally,  lead  to  the  conclusion  that  flow  modulation  has  no  effect  on  the 
average  behavior  of  the  flow  and  that  the  turbulent  structure  can  still  be 


described  by  quasi-steady  turbulence  models.  A  careful  study  of  these  in¬ 
vestigations,  however,  reveals  that  in  all  these  experiments,  the  frequency 
of  flow  modulation  was  small  compared  with  the  characteristic  frequency  of 
turbulence  in  the  boundary  layer.  One  can  expect  to  find  a  significant 
effect  of  flow  modulation  on  the  average  flow  structure  only  when  the  mod¬ 
ulation  frequency  is  comparable  to  the  characteristic  turbulent  frequency 
in  the  flow.  The  very  recent  report  of  Houdeville  and  Cousteix  (1978)  on 
the  unsteady  boundary  layer  in  the  neighborhood  of  separation  (where  this 
condition  is  satisfied)  does  indicate  strong  effects  of  flow  modulation  on 
the  flow  structure.  Hence,  the  behavior  of  an  unsteady  turbulent  flow  is 
not  only  determined  by  the  Strouhal  number  (ratio  of  the  time  scale  of  mean 
flow  to  the  time  scale  of  unsteadiness),  but  also  by  the  ratio  of  the  time 
scale  of  unsteadiness  to  the  characteristic  time  scale  of  turbulence.  Further, 
since  most  of  the  interaction  between  the  impressed  oscillation  and  the 
turbulence  structure  is  likely  to  occur  within  the  Stokes  layer,  it  is  neces¬ 
sary  to  design  the  experiment  carefully  so  as  to  get  a  thick  enough 
Stokes  layer  and  yet  keep  the  modulation  frequency  at  a  value  comparable  to 
the  characteristic  turbulent  frequency. 

The  present  study  was  directed  towards  obtaining  experimental 
data  on  fully  developed  oscillatory  pipe  flow  at  transitional  Reynolds  numbers. 
The  nature  of  the  periodic  flow  was  found  to  depend  strongly  on  the  nature 
of  the  steady  flow  at  the  mean  Reynolds  number.  It  was  found  that  if  the 
steady  mean  flow  was  fully  turbulent  (with  an  intermittency  of  unity),  the 
flow  would  remain  fully  turbulent  when  the  Reynolds  number  was  modulated 
at  a  frequency  large  enough  to  approach  the  characteristic  turbulent  frequency. 
On  the  other  hand,  if  the  steady  mean  flow  was  not  turbulent  but  contained 
turbulent  puffs  at  a  low  intermittency,  the  flow  would  tend  to  laminarize 
on  periodic  oscillation.  Both  these  flow  situations  were  studied  in  detail 
using  Laser  Doppler  anemometry  (LDA)  together  with  a  digital  phase  averaging 
technique.  Important  features  of  the  present  experiments,  particularly  with 
regard  to  the  turbulent  unsteady  flow  are  (i)  the  oscillation  frequency  was 
comparable  to  the  characteristic  turbulence  frequency  in  the  flow  and  (ii) 
the  viscous  sublayer  (2  mm)  and  Stokes  sublayer  (6  mm)  were  thick  enough  to 
permit  detailed  measurements  in  a  region  where  the  interaction  of  flow  oscil¬ 
lation  and  turbulence  have  been  found  to  be  very  strong. 


tet. 
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II.  EXPERIMENTAL  PROGRAM 

A.  General.  A  test  facility  for  producing  sinusoidally  oscil¬ 
lating  pipe  flow  at  transitional  Reynolds  numbers  was  designed  and  built 
for  the  present  study.  This  facility  used  Eureka  oil  (of  kinematic  vis¬ 
cosity,  v  =  1.384x10  5  m2/s  at  26°C)  and  consisted  of  a  long  pipe  of  50 

mm  diameter  in  which  the  oscillating  flow  was  maintained.  The  flow  oscil¬ 
lation  was  produced  by  varying  the  flow  exit  area  sinusoidally  by  a  rota¬ 
ting  sleeve  situated  at  the  exit  end  of  the  pipe.  The  design  of  the  rota¬ 
ting  sleeve  was  very  similar  to  that  used  in  the  classical  experiments  of 
Schultz-Grunow  (1940) .  The  exit  area  was  kept  small  in  comparison  with  the 
pipe  area  so  that  most  of  the  head  loss  in  the  system  occurred  at  the 
exit.  Under  such  conditions,  it  is  possible  to  obtain  a  discharge  rate  nearly 
proportional  to  the  exit  area.  Velocity  measurements  were  made  using  a 
single  channel  LDA.  The  data  was  acquired  and  processed  using  the  IIHR 
IBM/1800  Data  Acquisition  System  and  a  data  processing  technique  developed 
by  Mueller  and  Ramaprian  (1977)  for  the  analysis  of  unsteady  turbulent 
flows.  The  details  are  described  in  the  following  sections 

B,  Description  of  the  Apparatus.  The  experimental  apparatus  is 
shown  schematically  in  figure  1.  Oil  is  pumped  by  the  circulating  pump  from 
the  sump  tank  to  the  constant  level  overhead  tank.  The  overhead  tank  is  a 
cylindrical  reservoir  of  0.5  m  diameter  x  0.75  m  height.  The  head  is  main¬ 
tained  constant  (about  4  m  above  the  ground  level)  by  providing  a  50  mm 
diameter  overflow  pipe  which  returns  the  extra  oil  to  the  sump  tank.  Oil 
from  the  constant  head  tank  flows  down  through  the  100  mm  diameter  inlet  pipe 
and  then  through  the  200  mm  diameter  curved  pipe.  A  set  of  turning  vanes 

is  provided  in  the  curved  pipe  to  distribute  the  flow  uniformly  in  the  pipe. 
The  horizontal  section  of  the  200  mm  pipe  acts  as  a  calming  chamber.  This 
section  is  followed  by  a  bell-shaped  contraction  nozzle  of  50  mm  exit  dia¬ 
meter.  The  test  pipe  is  made  of  copper  and  has  a  nominal  inlet  diameter  of 
50  mm.  It  has  a  total  length  of  8.8  m,  (i.e.,  L/D  =  176)  which  ensures 
full  development  of  flow  even  under  laminar  flow  conditions.  The  pipe  has 
been  assembled  from  sections.  The  joints  between  the  different  sections 
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have  been  matched  with  special  care  insuring  that  there  are  no  irregulari¬ 
ties  in  the  surface.  Twenty-eight  wall  static  pressure  taps  of  1  mm  diame¬ 
ter  have  been  provided  at  regular  intervals  along  the  pipe  for  the  measure¬ 
ment  of  pressure  loss  through  pipe.  These  pressure  taps  have  been  provided 
along  the  bottom  of  the  tube  in  order  to  prevent  air  bubbles  from  accumula¬ 
ting  near  the  holes.  Provision  has  also  been  made  for  flush  mounting  pres¬ 
sure  transducers  at  three  locations  along  the  pipe.  The  first  one  is  at  a 
point  1.3  m  from  the  contraction  exit.  The  second  one  is  at  a  point  30  mm 
upstream  from  the  measurement  location  and  the  third  at  1.3  m  upstream  of 
the  second.  Pressure  transducers  mounted  at  locations  2  and  3  can  be  used 
to  obtain  the  pressure  gradient  in  the  pipe.  The  test  section  where  velo¬ 
city  measurements  were  made  is  a  plexiglass  tube,  0.3  m  long  and  50  mm  in 
internal  diameter.  The  upstream  end  of  the  tube  is  machined  and  matched 
with  the  copper  pipe.  The  other  end  is  connected  to  a  50  mm  diameter  x  300 
mm  long  copper  tube  whose  downstream  end  is  closed.  Two  longitudinal  rec¬ 
tangular  slots,  50  mm  x  3.2  mm,  milled  on  the  surface  of  this  tube  along  two 
generators  180  degrees  apart  from  each  other,  provide  the  exit  for  the  oil. 
The  oil  falls  into  a  rectangular  plexiglass  collection  chamber  (vented  ade¬ 
quately  to  the  atmosphere)  from  where  it  flows  into  a  weigh  tank.  The 
weigh  tank  rests  on  a  platform  weighing  machine  and  is  used  for  measuring 
the  mass  rate  of  flow  of  the  oil.  The  oil  from  the  weigh  tank  can  be 
drained  via  a  valve  into  the  sump  tank. 

The  technique  used  for  producing  discharge  oscillation  is  similar 
to  that  used  by  Schultz-Grunow.  The  device  consists  of  a  hollow  copper  cyl¬ 
indrical  sleeve  rotating  over  the  slotted  copper  tube  mentioned  above.  One 
end  of  the  sleeve  is  profiled  on  a  milling  machine  in  such  a  way  as  to  form, 
on  development,  two  cycles  of  a  sine  wave  of  amplitude  12  mm.  The  other  end 
of  the  sleeve  is  closed  and  is  attached  to  a  13  mm  diameter  shaft.  The 
shaft  is  connected  to  a  variable  speed  motor,  by  means  of  a  coupling  that 
permits  adjustment  of  the  longitudinal  position  of  the  sleeve  relative  to 
the  slots.  The  shaft  is  supported  by  a  bearing  housed  in  the  wall  of  the 
plexiglass  collection  chamber.  The  motor  is  mounted  independently  on  a 
bracket  secured  to  the  framework  supporting  the  apparatus.  The  motor  is  man¬ 
ufactured  by  Bodine  Electric  and  is  a  geared  D.C.  motor  having  an  output 
speed  range  of  1.2-55  rpm.  Its  speed  is  regulated  electronically  by 
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controlling  the  motor  voltage  via  a  feed-back  circuit.  The  regulator 
maintains  the  speed  of  the  motor  within  1%  of  the  preset  speed,  which 
is  continuously  selectable  within  the  range  mentioned  above.  When  the 
sleeve  rotates,  the  area  of  opening  of  the  pair  of  slots  is  continuously 
varied  by  the  profiled  sleeve.  It  is  seen  that  the  exit  area  variation  will 
be  a  sinusoidal  function  of  time  with  the  variation  going  through  two  cycles 
during  each  revolution  of  the  sleeve.  The  mean  exit  area  can  be  adjusted 
within  limits  by  changing  the  axial  location  of  the  sleeve  relative  to  the 
slots.  In  the  present  study,  this  location  has  been  chosen  (after  several 
trial  runs)  such  that  the  mean  flow  Reynolds  number  would  be  about  2100. 

The  sleeve  profile  was  initially  designed  for  an  estimated  discharge  modu¬ 
lation  amplitude  of  0.35  x  mean  discharge.  Later  tests  showed  that  the  ac¬ 
tual  discharge  modulation  at  infinitely  low  frequencies  (i.e.,  quasi-steady 
state),  was  in  fact,  quite  close  to  this  calue.  With  the  present  design, 
it  is  possible  to  produce  fully  laminar  flow  at  the  minimum  slot  opening  and 
fully  turbulent  flow  at  the  maximum  slot  opening  (under  steady  flow  conditions). 
At  higher  frequencies,  the  amplitude  is  reduced.  This  aspect  is  discussed 
next  in  greater  detail. 


C.  Performance  of  the  Apparatus 

1.  Quasi-steady  performance.  If  H  is  the  gross  head  available 
(about  2.7  m) ,  h^  the  overall  head  loss  in  the  system,  Vg  the  exit  velocity, 
V  the  test  section  velocity,  Ag  the  exit  area,  A  the  cross  sectional  area 
of  the  test  section  and  the  discharge  coefficient  (assumed  constant),  we 
have,  for  quasi-steady  flow,  2 

H  -  h  =  V  2/2g  -  V2/2g  =  (Since  V  is  very  much  (1) 

6  ®  smaller  than  V  ) 

e 


V  =  CD  -f  /2g(H  -  hf)  (2) 

If  h^  =  K  V2/2g,  where  K  is  a  constant  for  the  system,  one  finally  gets 

V  =  /2lH  (CDAe/A)  [1/U  +  ^(CpA^/A) 2  Y^] 

=  /2gif  CD  (Ae/A)F 


(3) 
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where 

F  =  [1/{1  +  K(C_A  /A)2}]15 
U  e 

Since  is  assumed  constant,  V  is  proportional  to  Ag,  if  F  is  constant; 
i.e.,  if  K(Ae/A)‘<<l.  The  frictional  loss  in  the  various  components  of 
the  system  can  be  estimated  using  standard  handbook  data.  These  estimates 
yield  a  value  of  Ks2  for  the  present  system.  The  ratio  Ag/A  has  a  maximum 
value  of  0.1  in  the  present  case  (for  fully  "open"  position).  Substituting 
these  values,  one  gets,  for  the  fully  "open"  sleeve  position, 

K(A  /A)2  =  0.02  (4) 

e 

This  means  that  F  varies  from  1.0  to  0.98  from  the  fully  closed  to 
fully  open  sleeve  positions.  Thus  an  almost  sinusoidal  discharge  rate  (or 
velocity)  can  be  expected.  The  quasi-steady  discharge  variation  with  the 
rotation  of  the  sleeve  is  shown  in  figure  2.  It  is  seen  that  the  modulation 
is  very  nearly  sinusoidal.  It  is  possible  to  improve  the  wave  shape  of 
modulation  by  departing  slightly  from  an  exact  sinusoidal  sleeve  profile. 

The  required  sleeve  profile  can  then,  be  calculated  from  Eq.  (3).  This 
requires  a  knowledge  of  the  precise  values  of  K  and  as  functions  of 
quasi-steady  discharge.  These  can  be  obtained  after  the  steady  flow  opera¬ 
tion.  In  the  present  design,  however,  this  procedure  was  not  followed. 

2.  Unsteady  performance.  The  performance  of  the  system  will 
be  described  by  Eq.  (3)  only  when  the  exit  area  is  constant  or  varied  infi¬ 
nitely  slowly  (i.e.,  under  quasi-steady  conditions).  At  finite  frequencies 
of  area  variation,  the  flow  in  the  system  will  be  governed  by  the  appropriate 
unsteady  euations  of  motion.  Specifically,  the  dynamic  response  of  the  system 
will  be  affected  by  the  inertia  of  the  fluid.  As  a  consequence,  as  the  fre¬ 
quency  of  oscillation  is  increased  the  amplitude  of  the  flow  modulation  will 
decrease,  the  modulation  will  be  distorted  from  the  sinusoidal  shape  and  the 
pressure  at  any  point  in  the  system  will  oscillate.  Further,  the  discharge 
and  pressure  oscillations  will  generally  be  out  of  phase  with  the  exit 
area  variation.  Some  of  these  effects  are  seen  from  figure  2  showing 
the  shape  of  the  discharge  (or  cross  sectional  average  velocity,  U^)  Vs.  time 
curve,  at  two  different  oscillation  frequencies;  0.057  Hz  and  1.75  Hz.  In 
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this  figure,  the  distribution  of  U  is  normalized  with  respect  to  the 

m 

quasi-steady  cross  sectional  average  velocity  at  the  mid  slot  opening. 

The  quasi-steady  velocity  distribution  (U  )  as  well  as  an  exact  sinu¬ 
soidal  distribution  are  also  shown  in  the  figure.  The  reduction  in  ampli¬ 
tude  and  the  phase  shift  of  discharge  modulation,  as  well  as  the  distortion 
in  the  shape  of  the  modulation  curve  can  be  easily  seen  from  the  figure. 

It  is  seen  that  at  the  lower  frequency,  these  effects  are  very  small  while 
at  the  higher  frequency,  they  are  perceptible.  We,  however,  regard  the  dis¬ 
tortion  from  the  sinusoidal  shape  to  be  still  small  enough  to  be  ignored 
for  the  purpose  of  the  present  study.  Likewise,  the  reduction  in  amplitude 
at  the  higher  frequency  is  no  cause  for  concern,  since  it  is  still  large 
enough  (20%  of  the  mean)  to  allow  useful  and  meaningful  measurements  to  be 
made  at  this  frequency. 

D.  Instrumentation 

1.  Pressure  measurement.  Wall  static  pressure  distribution  along 
the  pipe  in  steady  flow  was  measured  using  the  static  pressure  taps  and  an 
inverted  U-tube  oil  manometer.  The  most  upstream  static  pressure  hole  was 
connected  to  one  of  the  limbs  of  the  manometer.  The  other  limb  could  be 
connected  to  any  static  pressure  hole  via  a  manifold  with  which  all  the 
static  pressure  holes  were  in  communication  through  individual  stop  cocks. 

The  manometer  thus  measured  the  pressure  drop  from  the  first  hole. 

In  the  unsteady  flow  experiments,  pressure  was  measured  using  a 
pressure  transducer  mounted  flush  with  the  inner  surface  of  the  pipe.  Pres¬ 
sure  measurements  were  made  at  only  one  position  viz.  at  station  2  in  Fig.  1. 
Unsteady  pressure  drop  measurements  were,  however,  not  made  due  to  instrumen¬ 
tation  problems.  The  pressure  transducer  used  was  an  ENDEVCO  8510  semi¬ 
conductor  type,  of  range  0-15  psi,  with  a  sensitivity  of  9.5  mv/psi  at  an 
excitation  voltage  of  5  volts.  An  IIHR-built  stralngage  balance  and  5- 
volt-power  supply  was  used  with  the  transducer.  The  unbalanced  voltage 
output  from  the  bridge  was  amplified  by  a  PRESTON  D.C.  amplifier.  The 
instantaneous  amplified  output  which  was  proportional  to  the  static  pressure 
was  measured  using  the  data  acquisition  technique  described  later. 


2.  Discharge  measurements.  The  time  averaged  flow  rate  was 
measured  using  the  weigh  tqnk.  The  flow  rate  was  measured  by  clocking 
Che  time  required  for  a  preset  weight  of  oil  to  flow  into  the  weigh  tank. 
Unfortunately,  due  to  the  small  capacity  of  the  tank,  it  was  necessary  to 
limit  the  duration  of  oil  collection  to  15-20  s  only.  This  resulted  in  some 
inaccuracy  in  the  estimate  of  the  flow  rate,  especially  in  the  unsteady  flow 
experiments.  These  measurements  of  flow  rate  were,  however,  neither  cri¬ 
tical  nor  extremely  important  in  the  present  study.  In  fact,  the  flow  rate 
could  be  more  accurately  determined  in  specific  cases  by  integrating  the 
velocity  distribution  across  the  pipe  obtained  from  LDA  measurements. 

3.  Velocity  measurements.  Velocity  measurements  in  unsteady 
flow  studies  reported  so  far  in  the  literature  have  been  made  using  either 
pitot-tubes  or  hot-wire/film  anemometry.  These  probes  have  two  important 
disadvantages;  (i)  they  disturb  the  flow  and  (ii)  they  do  not  allow  reversals 
in  velocity  to  be  detected  or  measured.  In  addition,  the  Pitot  tube  has 

the  added  limitations  that  it  cannot  measure  turbulent  intensity  nor  can 
it  even  by  used  in  any  except  very  slowly  oscillating  flows.  One  of  the 
most  important  features  of  the  present  study  is  the  use  of  LDA,  a  technique 
that  can  measure  instantaneous  flow  velocities  without  disturbing  the 
flow.  The  LDA  has  also  the  additional  advantages  that  its  signal  is  propor¬ 
tional  to  the  velocity  and  that  it  does  not  require  calibration  against 
other  standards.  One  can,  thus,  avoid  the  calibration  drift  problems 
associated  with  hot-film  anemometry  in  liquids.  Reversal  of  flow  can  be 
detected  and  measured  using  LDA  if  a  "frequency  shift"  technique  is  used. 

The  LDA  system  used  in  the  present  studies  used  a  frequency-shift  feature 
though  flow  reversal  did  not  actually  occur  in  the  range  of  flows  studied. 

The  single  channel  LDA  system  of  11HR  and  the  details  of  the  traverse 
system  are  described  in  appendix  A. 

A.  Reference  time  signal.  In  the  present  report,  the  time  his¬ 
tory  of  the  flow  during  the  oscillation  cycle  was  studied  using  the  instant 
of  maximum  slot  opening  as  the  reference  (t  =  0)  point  in  the  cycle.  As 
will  be  seen  later,  the  Data  Acquisition  System  was  designed  to  use  this 
point  as  the  instant  at  which  data  sampling  should  begin.  For  this  purpose 
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an  electromagnetic  device  was  used  to  produce  a  voltage  pulse  at  the 
instant  of  maximum  slot  opening  once  during  each  revolution  of  the  pro¬ 
file  sleeve.  (Note  that  this  corresponds  to  once  in  every  two  oscillation 
cycles).  The  device  included  an  electromagnetic  pick  up  (1.5  nan)  diameter) 
mounted  on  the  support  frame  and  a  thin  steel  rod  (1.5  mm  diameter)  atta¬ 
ched  to  the  motor  shaft.  The  position  of  the  rod  was  carefully  adjusted 
such  that  during  its  rotation  the  tip  passed  within  1  mm  of  the  pick  up, 
cutting  the  magnetic  field.  The  voltage  induced  in  the  pick  up  was  used 
as  the  reference  time  signal.  The  intensity  and  duration  of  the  pulse 
would  depend  upon  the  diameter  of  the  rod  and  its  rotational  speed.  The 
diameter  of  the  rod  (1.5  mm)  was  determined  from  the  requirement  that  a  sig¬ 
nal  of  sufficient  strength  to  drive  the  Data  Acquisition  System  was  to  be 
produced  at  the  lowest  speed  studied.  With  this  diameter  and  with  the  radius 
of  rotation  of  the  rod  of  about  150  mm,  the  pulse  duration  was  found  to 
vary  from  about  2.5  degrees,  at  the  lowest  speed,  to  about  3.5  degrees  at 
the  highest  speed  of  rotation.  This  corresponds  to  the  range  5-7  degrees 
when  referred  to  an  oscillation  cycle  and  represents  the  overall  uncertainty 
in  the  reference  time  signal.  The  phase  angle  measurements  reported  in  this 
study  are  all  subjected  to  this  uncertainty.  The  angular  position  of  the 
rod  was  adjusted  such  that  the  rod  was  aligned  with  the  pick  up  in  the 
maximum  slot  opening  position.  A  360-degree  graduated  protractor  fixed  to 
the  plexiglass  collection  chamber  and  a  pointer  attached  to  the  shaft 
enabled  setting  up  the  sleeve  at  any  desired  position  (6)  relative  to 
the  maximum  slot  opening  position  (t  =  0  or  9  =  0)  in  steady  flow  experi¬ 
ments. 


E.  Experimental  Details 

1.  General.  As  already  mentioned,  two  series  of  experiments 
were  carried  out.  In  both  the  series,  the  steady  flow  had  a  Reynolds 
number  of  about  2900  at  maximum  slot  opening  (0=0°)  and  a  Reynolds  number 
of  about  1300  at  minimum  slot  opening  (0=180°).  The  mean  steady  flow  Rey¬ 
nolds  number  was  about  2100  and  this  occured  at  0s9O°.  However,  in  the 
first  series,  the  steady  flow  at  0=90°  was  found  to  be  fully  turbulent  at 
all  times.  In  fact,  the  intermittency  of  turbulence  remained  at  unity 
for  6S  100°.  When  this  flow  was  oscillated  at  the  highest  possible  frequency 


of  1.75  Hz,  the  flow  remained  fully  turbulent.  Instantaneous  velocity 
measurements  were  obtained  across  the  pipe  for  this  situation  (Run  13).  A 
rough  estimate  of  the  turbulent  burst  frequency  in  this  flow  using  the  cri¬ 
terion  of  Rao  et  al,  (1971)  indicates  a  value  of  about  2  Hz.  Thus  the 
oscillation  frequency  can  be  expected  to  interact  significantly  with  the 
turbulent  structure  in  this  flow.  In  the  second  series  of  experiments,  while 
the  steady  flow  at  the  maximum  and  minimum  slot  opening  behaved  exactly  as 
in  the  first  series,  the  mean  flow  at  0=90°  exhibited  an  intermittent  turbu¬ 
lent  structure.  The  structure  strongly  resembled  the  puff  type  transitional 
structure  studied  by  Wygnanski  and  Champagne  (1973).  In  fact,  the  steady 
flow  became  fully  turbulent  only  at  values  of  0<6O°.  When  this  flow  was 
oscillated  (at  whatever  frequency)  it  was  found  to  laminarize  with  the  inter- 
mittency  of  puffs  dropping  almost  to  zero.  This  is  clearly  seen  from  figure  3 
where  a  photograph  of  the  storage  oscilloscope  traces  corresponding  to  the 
center  line  velocity  signals  obtained  from  the  LDA  under  three  different 
conditions  are  shown  superimposed  on  one  another.  These  are:  i)  steady 
flow  at  0=0°,  Re  =  2872  (fully  turbulent),  ii)  steady  flow  at  6=180°,  Re  = 

1278  (fully  laminar)  and  (iii)  oscillatory  flow  at  0.057  Hz  between  the 
above  two  Reynolds  numbers  (laminarized) .  The  peak  velocity  in  the  os¬ 
cillatory  flow  is  different  from  the  velocity  in  steady  turbulent  flow 
because  the  steady  flow  is  laminar  and  hence  has  a  different  velocity 
distribution.  Velocity  measurements  were  made  in  this  laminarized 
unsteady  flow  at  oscillation  frequencies  of  1.75  Hz  (Run  23)  and  0.057 
Hz  (Run  24).  These  frequencies  correspond  to  Strouhal  numbers  S  (defined 
as  --S-P-)  of  1.0  and  0.032  respectively 

In  both  the  series,  measurements  of  velocity  distribution  in  the 
pipe  were  made  for  the  steady  flows  at  0=0  and  0=180°.  In  addition,  steady 
flow  axial  pressure  drop  data  were  obtained  for  several  (fixed)  angular 
positions,  0  of  the  sleeve.  These  steady  flow  experiments  can  be  consid¬ 
ered  to  approximate  closely  to  quasi-steady  experiments.  They  are  used 
as  a  reference  in  the  study  of  the  unsteady  flow  behavior.  The  use  of 
these  data  instead  of  steady  flow  data  from  other  sources  is  both  necessary 
and  desirable  in  view  of  the  very  low  Reynolds  numbers  and  transitional 
nature  of  the  flow. 
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2.  Quasi-steady  measurements.  One  of  the  features  of  tne 
present  design  is  that  it  allows  one  to  make  measurements  at  steady  flow  con¬ 
ditions  corresponding  to  any  desired  point  on  the  oscillation  cycle.  These 
measurements  closely  approximate  to  quasi-steady  flow  measurements,  i.e., 
measurements  obtained  in  unsteady  flow  oscillating  at  infinitesimally  low 
frequencies.  These  quasi-steady  measurements  form  the  basis  for  studying 
the  effect  of  flow  oscillation  at  finite  frequencies.  Quasi-steady 
measurements  were  made  in  both  series  of  experiments.  These  included: 

(1)  discharge 

(2)  pressure  drop  and 

(3)  velocity  distribution  across  the  pipe. 

The  discharge  measurements  were  made  using  the  weigh  tank  and 
have  been  already  referred  to  earlier  in  Section  D.2.  As  mentioned  pre¬ 
viously,  it  was  found  that  due  to  the  small  collection  time  used,  there  was 
some  inaccuracy  introduced  into  the  discharge  measurements.  The  discharge 
data  obtained  from  the  weigh  tank  measurements  were  compared  with  the 
more  accurate  values  obtained  from  the  integration  of  velocity  profiles 
in  the  test  section  and  a  correction  factor  for  the  former  was  arrived  at. 
This  was  used  for  the  rest  of  the  discharge  data  from  the  weigh  tank.  These 
discharge  measurements  were  made  for  several  fixed  slot  openings  (i.e.,  for 
several  fixed  values  of  0)  in  both  the  experimental  series. 

The  distribution  of  static  pressure  along  the  pipe  was  measured 
only  in  series  2.  Pressure  drop  measurements  were  made  for  several 
values  of  9.  These  pressure  measurements  were  made  using  the  inverted 
oil  manometer  mentioned  earlier. 

The  velocity  measurements  were  made  in  both  series.  However, 
they  were  made  only  for  two  values  of  9;  9=90°  and  0=0°  (Runs  11,  12;  21, 

22).  Further,  the  velocity  profiles  were  measured  at  16  points  across  the 
pipe  in  series  1.  In  series  2,  velocities  were  measured  at  a  fewer  number 
of  points.  Since  these  measurements  appeared  to  reproduce  closely  the 
first  series  of  measurements,  it  was  felt  unnecessary  to  obtain  data  at 
finer  intervals.  The  velocity  measurements  were  made  using  the  IIHR 
Data  Acquisition  System,  as  described  in  Section  F.  Measurements  at  9=90° 
(i.e.,  in  laminar  flow)  served  to  check  the  accuracy  of  measurement  as 
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well  as  the  quality  of  flow  (axisymmetry)  in  the  pipe  by  comparing  the  ob¬ 
servation  with  the  theoretical  parabolic  profile.  Measurements  in  fully 
turbulent  flow  (0=0°)  provided  the  basis  for  studying  the  effect  of  flow 
oscillation.  The  turbulent  velocity  distribution  in  the  pipe  at  low  Rey¬ 
nolds  numbers  is  strongly  dependent  on  the  Reynolds  number  and  possibly  on 
the  characteristics  of  the  specific  experimental  apparatus.  It  is  thus 
necessary  to  use  the  quasi-steady  measurements  in  the  same  apparatus  at 
similar  Reynolds  numbers  in  order  to  make  any  comparisons  with  unsteady 
flow  measurements. 

3.  Unsteady  flow  measurements.  The  flow  was  oscillated  at  the 
desired  frequency  by  setting  the  voltage  control  of  the  motor  at  the  appro¬ 
priate  position.  The  rotational  speed  (equal  to  twice  the  oscillation  fre¬ 
quency)  was  measured  by  clocking  the  time  taken  for  50  revolutions  of  the 
sleeve.  Several  minutes  were  allowed  to  elapse  after  starting  the  oscilla¬ 
tions  before  the  measurements  were  made.  As  in  the  case  of  quasi-steady 
flow,  velocity  measurements  were  made  at  16  points  across  the  pipe.  The 
instantaneous  velocity  signal  from  the  LDA,  was  read  on  to  he  IIHR  Data 
Acquisition  System  and  processed  using  the  procedure  described  in  section  F. 
The  reference  time  signal  from  the  magnetic  pick-up  was  used  as  the  syn¬ 
chronizing  signal  in  the  data  acquisition  process. 

An  oscillation  frequency  of  1.75  Hz  was  studied  in  series  1. 

(Run  13).  As  already  mentioned,  the  flow  was  observed  to  be  fully  turbu¬ 
lent  in  this  case.  At  lower  frequencies  (say  <1  Hz),  the  flow  was  found  to 
switch  intermittently  between  laminar  and  turbulent  regimes.  This  was 
accompanied  by  large  variations  in  the  "mean"  velocity  of  the  flow.  Hence, 
it  was  not  possible  to  make  any  quantitative  velocity  measurements  under 
these  conditions.  In  the  second  series  of  experiments,  the  flow  remained 
laminar  or  nearly  laminar  at  1.75  Hz  (Run  23)  as  well  as  at  a  very  low 
frequency,  viz.  0.057  Hz  (Run  24).  Detailed  velocity  profiles  were  obtained  at 
these  two  frequencies. 

4.  Intermittency  measurements.  A  careful  and  extended  study  did 
reveal  that  even  in  the  second  series  of  experiments,  the  flow  did  not  remain 
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fully  laminar  at  all  oscillation  frequencies.  It  became  turbulent  inter¬ 
mittently  and  this  intermittency  was  found  to  depend  on  the  oscillation 
frequency.  In  order  to  study  this  aspect  in  some  detail,  the  velocity  sig¬ 
nal  from  the  LDA  was  recorded  on  a  strip  chart  recorder  for  a  duration  long 
enough  to  allow  one  to  estimate  the  intermittency  of  turbulence  from  the 
record.  This  duration  corresponded  approximately  to  300  cycles  or  5  minutes 
whichever  was  larger.  Such  records  were  obtained  for  several  oscillation 
frequencies  in  the  range  of  0.05  Hz  -  1.75  Hz  and  for  two  points  in  the  pipe, 
viz . ,  r  =  0  mm,  and  r  =  18. 4  mm. 


F.  Data  Acquisition  and  Processing.  The  procedure  used  in  the 
present  study  for  the  analysis  of  unsteady  turbulent  flows  is  based  on  that 
used  by  Nakato  (1974)  and  is  essentially  similar  to  that  used  by  many  recent 
investigators  of  oscillatory  turbulent  boundary  layers.  The  author  consid¬ 
ers  it  to  be  superior  and  more  informative  than  the  analog  procedures  used 
by  previous  investigators  of  pulsating  flows  in  arteries  or  rigid  walled- 
tubes.  The  analysis  can  be  used  in  a  degenerate  form  for  unsteady  laminar 
flows,  steady  turbulent  and  steady  laminar  flows  also.  The  analysis  is 
briefly  described  below. 

Let  the  instantaneous  velocity  U (y ,  ,t)  be  written  as: 

U(r,0,t)  =  U(r)  +  U  (r,0)  +  u(r,6,t)  =  <  U(r,0)  >  +  u(r,6,t)  (8) 

P 

where  U  is  the  time  mean  velocity,  u  is  the  turbulent  (random)  velocity, 
and  <U>  is  the  deterministic  velocity.  The  deterministic  velocity  can  be 
obtained  by  a  process  of  ensemble  averaging.  In  the  present  periodic  flow, 
ensemble  averaging  is  equivalent  to  phase  averaging  i.e.,  averaging  over 
the  values  obtained  at  identical  values  of  r  and  0  in  a  large  number  of 
cycles.  Thus, 


<  U (  r, 0)  >  =  -  ^  U.(  r,0,t)  (9) 

where  N  is  the  total  number  of  cycles  being  averaged  and  IT  is  the  velo¬ 
city  in  the  ith  cycle.  The  time  mean  value  U(r)  is  obtained  by  long  time 
averaging  or  equivalently  by  averaging  <U>  (r,u)  over  a  cycle.  Thus, 


U(r )  =  y-  /  <U  (  r,  0)>dG 


If  M  descrete  samples  are  taken  during  a  cycle,  we  can  write 


(10) 
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(11) 

(12) 

(13) 

(14) 

(15) 


u’<r)  =/H  jSi  <0^) 


In  the  present  experiments,  the  Data  Acquisition  System  was  programmed  to 
sample  the  LDA  output  once  at  the  end  of  each  sampling  interval  (AT), 
starting  from  the  instant  when  the  synchronizing  signal  (reference  time 
signal)  was  received.  AT  was  so  chosen  that  100  ATS  the  period  of  oscil¬ 
lation.  This  would  ensure  that  100  samples  could  be  taken  per  cycle.  How¬ 
ever,  only  96  samples  were  taken,  no  samples  being  taken  at  the  end  of 
97th,  98th,  99th,  and  100th  time  intervals.  These  "holes"  in  the  data 
were  required  in  order  to  allow  for  the  fact  that  the  period  of  oscilla¬ 
tion  might  not  be  exactly  equal  to  but  less  than  100  AT.  It  may  be  noted 
that  with  the  present  program,  data  would  be  acquired  actually  from  two  con¬ 
secutive  oscillation  cycles,  the  first  cycle  yielding  50  (No.  1  -  No.  50) 
samples  and  the  second  cycle  yielding  47  (No.  50  -  No.  96)  samples.  A 
total  of  100  revolutions  were  used  for  phase  averaging  (N  =  100)  at  the 
higher  oscillation  frequency.  At  the  lower  frequency,  25  revolutions  were 

used.  From  the  sampled  data,  the  various  velocities  U,  U  ,  u',  etc  were 

P 

computed  according  to  Eqs  (8)  through  (12).  In  most  cases,  the  experiments 
were  repeated  thrice  and  the  results  averaged  over  the  three  experiments. 

If  the  results  in  consecutive  experiments  differed  significantly  from  one 
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another,  the  experiments  were  repeated  till  three  consecutive  experiments 
gave  nearly  the  same  results.  (Such  repititions  were,  however,  rarely 
necessary) .  The  value  of  U  varied  insignificantly  among  consecutive 
repititions,  while  the  scatter  in  <U>  was  within  1%  and  the  scatter  in 
u'  within  5%.  The  computer  programs  in  Assembler  Language  used  in  the  pres¬ 
ent  experiments  for  both  the  acquisition  and  the  subsequent  processing  of  the 
data  were  developed  by  Mueller  (see  Mueller  and  Ramaprian  1977)  from  an 
earlier  version  developed  by  Nakato  (1974).  A  listing  of  the  program 
package  used  in  the  present  experiments  is  given  in  Appendix  C.  One  of 
the  important  features  of  the  present  program  was  that  it  would  detect 
the  LDA  signal  drop-out  and  reject  the  samples  taken  during  such  a  cycle. 

Data  acquisition  would  then  proceed  until  N  valid  revolutions  were  sampled. 
However,  with  the  use  of  frequency  shift  technique,  there  was  virtually  no 
"signal  dropout"  problem  and  this  feature  in  the  program  was  never  made  use  of. 

The  above  program  package  was  also  used  for  processing  unsteady 
laminar  flows,  the  value  of  u'  in  this  case  being  used  as  a  measure  of  the 
noise  in  the  optics  and  electronics.  It  was  found  that  u'  ranged  from 
2-5  parts  per  1000  in  these  cases,  indicating  the  acceptability  of  the 
experimental  procedure.  The  program  package  was  also  used  for  processing 
steady  flow  (laminar  and  turbulent)  data,  the  quantities  0  and  u'  being  the 
only  quantities  sought  for  in  these  cases.  While  experimenting  with  steady 
flows,  a  fictitious  synchronizing  signal  was  provided  from  a  square  wave 
generator  with  a  period  of  approximately  3  seconds.  With  100  "cycles"  of 
sampling,  this  would  result  in  an  overall  averaging  time  of  5  minutes. 

The  unsteady  data  acquisition  program  package  was  also  used  for 
the  unsteady  pressure  measurements  referred  to  in  section  D.l. 


III.  RESULTS  AND  DISCUSSION 

A.  Steady  Flow  Measurements 

1.  Axial  pressure  drop.  The  nondimens ional  pressure  drop 

AP*  defined  by 

AP*  =  (P . -P) / (1/2  pU  2)  (17) 

A  qm 


is  shown  in  figure  4  for  various  values  of  0  (i.e.,  for  various  Reynolds 
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numbers).  P  is  the  static  pressure  at  an  upstream  reference  point,  x=0 

A 

(see  figure  1).  P  is  the  static  pressure  at  a  longitudinal  location  x  and 
p  is  the  fluid  density.  The  values  shown  are  from  the  second  series  of 
experiments.  The  figure  shows  that  mean  flow  has  attained  full  development 
at  all  Reynolds  numbers  for  x*  =  x/D  >  70.  The  quasi-steady  friction 
factor,  X  is  defined  in  the  usual  way  as 

q 


rdP 


X  (9)  =  -  (6)  D] /[  1/2  pU 


r,  2i  = 


(18) 


q  '  '  *  dx . q  m  '  dx* 

and  thus,  is  given  by  the  slope  of  the  lines  in  figure  4.  The  values  of 
X^  at  different  Reynolds  numbers  are  shown  in  figure  5.  Also  shown  for 
comparison  are  the  laminar  relation,  X  =  64/Re  and  the  Blasius  relation 

_p. 

for  low  Reynolds  number  turbulent  flows,  X  =  0.3165  Re  u  (Schlichting  1968). 
The  value  of  used  in  figures  4  and  5  were  obtained  from  discharge 

measurements  made  using  the  measuring  tank.  Figure  5  is  presented  to 
show  that  the  apparatus  was  functioning  as  expected  and  that  the  flow  was 
transitional. 


2.  Velocity  distribution.  The  results  of  velocity  measurement 
at  0=180°  (minimum  slot  opening)  and  0=0°  (maximum  slot  opening)  are  shown 
in  figure  6.  In  each  case,  data  obtained  in  both  the  series  of  experiments 
are  shown  even  though  fewer  points  are  shown  for  series  2.  This  is  con¬ 
sidered  adequate  in  view  of  the  excellent  repeatability  observed  between 
the  two  series  of  steady  flow  experiments.  In  the  case  of  laminar  flow 
(0=180°),  the  data  are  seen  to  agree  well  with  the  theoretical  parabolic 
profile.  This  good  agreement,  however,  has  been  obtained,  as  already 
mentioned,  by  applying  a  constant  "zero  correction"  to  the  apparent  loca¬ 
tion  of  the  measurement  point.  The  adjustment  is  justified  by  the  fact  that 
a  constant  zero  correction  brings  the  measured  profile  in  agreement  with 
the  parabolic  profile.  The  same  correction  was  applied  to  the  turbulent  flow 
(0=0°)  measurements  also.  It  is  seen  from  the  figure  that  the  mean  velocity 
profiles  in  turbulent  flow  obtained  from  the  two  series  of  experiments  also 
indicate  agreement  with  each  other.  The  agreement  between  the  two  sets  of 
data  obtained  with  a  gap  of  several  months  in  between  establishes  the  accuracy 
of  the  measurement  procedure.  It  also  confirms  that  the  behavior  of  the 


18 


quasi-steady  flow  at  both  the  maximum  and  minimum  end  of  the  oscillation 
cycle  did  not  change  even  though  the  transitional  character  at  intermediate 
Reynolds  numbers  had  changed. 

The  velocity  profile  for  the  turbulent  flow  at  0=0°  is  shown  in 
figure  7  in  the  usual  wall  layer  coordinate  U+  ^U/u*)  vs.  Y+[=(D/2-r)u^/v] 

u  being  the  shear  velocity,  A  /p  where  x  is  the  wall  shear  stress.  The 

*  w  w 

value  of  u^  was  obtained  from  the  quasi-steady  pressure  drop  and  discharge 
data  for  0=0°  using  the  relation. 


X 

q 


(19) 


The  measured  velocity  distribution  is  compared  with  the  universal  law  of  the 
wall,  namely  U+  =  Y+  in  the  viscous  sub-layer  and  lT+  =  5.75  log  Y+  +  5.5 


in  the  fully  turbulent  region.  The  deviation  of  the  experimental  data  from 


the  universal  log  law  is  to  be  expected  in  view  of  the  very  low  Reynolds 


number  of  the  flow.  On  the  other  hand,  the  fact  that  the  first  data  point 
near  the  wall  falls  on  the  linear  U+  =  Y+  curve  indicates  not  only  that  this 


point  is  in  the  viscous  sub-layer  but  also  that  the  value  of  u^  and  r  are 
accurate. 


3.  Turbulence  intensity  distribution.  Figure  6  introduced 
earlier,  also  shows  the  distribution  of  the  rms  intensity,  u'  of  the 
longitudinal  turbulent  velocity  fluctuations  in  the  steady  turbulent  flow 
at  9=0.  The  distribution  is  normalized  using  the  velocity  and  length 
scales,  u^  and  D/2  respectively.  Consequently,  the  distribution  can  be 
expected  to  show  Reynolds  number  dependence  in  the  inner  region.  This  is, 
in  fact,  the  case  as  is  seen  from  the  comparison  with  the  data  of  Laufer 
(1954)  obtained  at  a  Reynolds  number  of  5xl05.  The  distribution  in 
the  region  0.5  <(l-2r/D)<l  seen  to  coincide  reasonably  well  with  Laufer' s 
data  while  in  the  region  (1—2 r /LJ><  0.5,  a  strong  Reynolds  number  effect  is 
observed.  In  fact,  the  large  viscous  region  in  the  present  case  allows 
one  to  observe  the  gradual  decrease  of  u'  from  a  maximum  towards  zero  at 
the  wall,  a  feature  which  is  usually  difficult  to  observe  in  higher  Rey¬ 
nolds  number  laboratory  flows.  The  purpose  of  presenting  the  steady  flow 
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turbulence  distribution  in  figure  6  is  to  provide  a  basis  for  the 
comparison  with  the  unsteady  turbulent  flow  data.  It  is  clear  that 
it  is  very  important  to  make  comparisons  in  approximately  the  same  Rey¬ 
nolds  number  range  in  order  to  arrive  at  reliable  conclusions. 


B.  Velocity  Distribution  in  Unsteady  Flow.  Figure  8  shows 
the  distribution  of  the  time-averaged  velocity  U  across  the  pipe  in  the 
three  runs  13,  23,  and  24.  The  velocities,  in  each  case,  are  normalized 
with  respect  to  the  corresponding  velocity  at  the  pipe  axis  (umax) •  The 
distributions  are  compared  with  the  quasi-steady  velocity  distributions 
corresponding  to  0=180°  and  6=0°.  It  is  seen  from  the  figure  that  the  un¬ 
steady  flows  in  run  23  and  run  24  that  appeared  to  be  laminar  were,  indeed, 
very  nearly  laminar  i.e.,  they  exhibit  nearly  parabolic  velocity  distribu¬ 
tions.  In  fully  developed  oscillatory  laminar  pipe  flow  of  a  Newtonian 
fluid,  the  time-mean  velocity  distribution  should  be  identical  to  that  of 
the  steady  Poiseuille  flow  since  the  velocity  field  is  determined 
from  the  linear  equation 

3U  1  3P  x  r  32U  1  3U. 

3F=  “7  9^+  V[a^+  r  3?]  (20) 

It  is  seen  that,  in  the  high  frequency  run  23,  in  which  the  inter- 
mittency  of  turbulent  puffs  was  negligibly  small  (see  figure  23),  the 
time-mean  unsteady  flow  behaves  almost  exactly  like  laminar  flow.  The 
deviation  of  the  experimental  data  of  run  24  from  the  parabolic  profile 
has  been  caused  by  the  flow  becoming  turbulent  with  a  larger  intermittency 
than  in  run  23.  It  is,  however,  important  to  observe  from  figure  8  that 
the  effect  of  flow  oscillation  on  the  time-mean  velocity  gradient  and  hence  the 
wall  shear  stress,  is  negligible  in  both  these  runs. 

When  the  oscillatory  flow  is  turbulent,  the  total  shear  stress 
(x)  distribution  in  the  pipe  is  given  by  the  linear  equation 

3U  _  1  3P  1  3(tr)  (21) 

3t  p  3x  pr  3r 

and  hence  the  distribution  of  time  averaged  shear  stress,  t  is  now  given 

by 

t/p  =  v§£  -  uv  (22) 


(where  u  and  v  are  the  turbulent  velocity  fluctuations  in  the  x  and  r 
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directions,  respectively),  the  velocity  field  is  not  determined  by  a 
linear  equation.  The  time  averaged  velocity  distribution,  and  hence  the 
time  averaged  wall  shear  stress  need  not,  therefore,  be  necessarily  the 
same  as  in  the  steady  flow  at  the  mean  Reynolds  number.  It  is,  in  fact, 
seen  from  figure  that  the  time-mean  velocity  distribution  in  the  unsteady 
turbulent  flow  of  run  13  exhibits  a  point  of  inflection  near  the  wall  and 
thus  differs  from  the  steady  state  turbulent  velocity  profile.  The  time 
mean  velocity  gradient  at  the  wall  is  slightly  larger  than  in  the  steady 
case  resulting  in  a  larger  time-mean  wall  shear  stress.  This  observation 
differs  from  those  reported  by  previous  investigators  of  oscillatory  tur¬ 
bulent  boundary  layers,  who  concluded  that  unsteadiness  has  no  effect  on  the 
time-mean  properties  (or  even  the  turbulent  structure)  of  the  flow.  How¬ 
ever,  as  already  mentioned,  a  significant  difference  between  the  present 
experiment  and  the  earlier  experiments  is  that  the  oscillation  frequency  in 
the  present  case  is  comparable  to  the  characteristic  turbulent  frequency. 

Figure  9  presents  the  variation  of  the  periodic  component  of 

the  velocity  through  the  oscillation  cycle  at  a  few  typical  points  across 

the  pipe.  The  velocity  is  normalized  with  respect  to  the  amplitude  (U  ) 

v  v  J  qmp  max 

of  the  quasi-steady  cross  sectional  average  periodic  velocity.  The  varia¬ 
tions  are  shown  only  for  700  degrees  of  oscillation  cycle,  short  record 
gaps  being  present  (for  reasons  mentioned  earlier)  from  0  -  7.2°  and  698.4 
-  720°,  i.e.,  a  total  of  4  sampling  periods.  For  run  13,  the  data  are  shown 
only  for  600  degrees  of  oscillation  cycle,  due  to  some  difficulties  encoun¬ 
tered  during  data  acquisition.  These  difficulties  were  subsequently  removed 
before  the  second  series  of  experiments  were  started.  However,  since  data 
over  only  360  degrees  of  oscillation  cycle  are  sufficient  for  the  purpose  of 
analysis,  the  loss  of  redundant  data  is  not  serious.  In  each  of  the  figures 
9  (a),  (b) ,  (c) ,  the  distributions  of  the  unsteady  cross  sectional  average 


periodic  velocity. 


U  and  quasi-steady  cross 


sectional  average  periodic 


velocity,  U  (proportional  to  the  effective  modulation  in  exit 

qmp 


area)  are  also  shown  for  comparison.  The  latter  curve  is  the  same  as  the 


periodic  part  of  the  velocity  distribution  shown  in  figure  2.  A  comparison 


of  the  laminar  flow  data  for  runs  23  and  24  shows  that  the  phase-wise  varia¬ 


tions  are  strongly  dependent  on  the  oscillation  frequency.  The  variations 


21 


appear  to  be  very  nearly  sinusoidal  at  the  lower  frequency  while  signi¬ 
ficant  distortion  in  the  wave  shape  can  be  observed  at  the  higher  oscil¬ 
lation  frequency.  The  distortion  at  the  higher  frequency  is,  as  mentioned 
earlier,  due  to  nonlinear  effects. 

It  is  also  clear  from  the  figures  9  (a) ,  and  9  (b) ,  that  at  the 
higher  oscillation  frequency,  there  is  a  significant  phase  shift  between 
the  velocity  variation  and  the  exit  area  opening.  Most  of  this  phase  shift 
is  brought  about  by  the  global  dynamics  of  the  system  and  is  observable  as 
the  phase  shift  of  unsteady  cross  sectional  average  flow.  However,  there 
is  also  a  relative  phase  shift  between  the  local  velocity  and  the  cross 
sectional  average  velocity.  The  flow  in  the  wall  region  leads  the  cross 
sectional  average  flow  while  the  flow  in  the  core  region  lags  behind  the 
average  flow.  Unfortunately,  owing  to  the  uncertainty  in  phase  angle  measure¬ 
ment  (±  3.5°)  and  the  distortion  in  the  modulation  profile,  it  is  not 
possible  to  make  a  precise  quantification  of  the  phase  differences  nor  is 
it  possible  to  look  for  differences  between  the  laminar  and  turbulent  flow. 

Figures  10  (a),  (b) ,  and  (c)  show  cross  plots  of  the  distribution 
of  the  periodic  velocity,  across  the  pipe  for  a  few  typical  fixed  phase 
angles  in  the  oscillation  cycle.  It  is  seen  that  at  the  higher  oscillation 
frequency,  the  variations  in  are  confined  essentially  to  q  =  (l-2r/D) < 0. 25 
(i.e.,  the  Stokes  layer),  the  rest  of  the  flow  oscillating  virtually  as 
one  solid  mass  at  all  times.  At  the  lower  frequency,  U  varies  across  the 
entire  pipe  at  all  times.  The  responses  of  turbulent  and  laminar  flows 
appear  to  be  qualitatively  similar. 

Figure  11  shows  the  distribution  of  the  "amplitude",  (U  )  of 
°  p  max 

the  periodic  velocity  across  the  pipe  for  the  three  experiments.  Since 
the  velocity  modulation  is  distorted  from  a  sine  wave  especially  at  the 
higher  frequency,  it  is  not  strictly  appropriate  to  use  the  term  "ampli¬ 
tude"  without  making  a  formal  harmonic  analysis.  However,  in  the  present 

case,  (U  )  denotes  one  half  of  the  peak  to  peak  variation  in  the  periodic 
p  max 

velocity  Up.  In  the  figure  (Up)max  is  normalized  with  respect  to  the  ampli¬ 
tude  (U  )  of  the  cross  sectional  average  periodic  velocity.  It  is  more 
mp  max 

clearly  seen  from  this  figure  that  at  the  higher  frequency,  the  amplitude 
variations  are  confined  to  a  small  region (n  <  0.25)  near  the  wall  while  at 
the  lower  frequency  the  variations  extend  over  the  entire  pipe  cross  section. 
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Also,  it  is  seen  that  an  overshoot  occurs  in  the  amplitude  distribution. 

The  dominant  effect  of  oscillation  frequency  is  clearly  seen  from  this 
figure.  At  the  same  time,  it  is  possible  to  notice,  in  this  figure,  a  dif¬ 
ference  between  the  responses  of  laminar  and  turbulent  flow  in  the  near¬ 
wall  region.  The  turbulent  flow  data  of  run  13  exhibit  a  larger  overshoot 
than  the  laminar  flow  data  of  run  23.  Also  shown  in  the  figure  are  the  results 
computed  from  the  exact  solution  of  Uchida  (1956)  for  fully  developed 
laminar  periodic  flow  in  a  pipe  under  the  influence  of  a  pressure  gradient 
sinusoidally  varying  with  time.  The  pressure  gradient  variation  in  the 
present  experiments  at  the  higher  frequency  differs  significantly  from  that 
assumed  in  Uchida' s  analysis.  The  experiment  at  the  lower  frequency,  how¬ 
ever,  approximates  closely  to  the  theoretical  case.  The  data  for  this  ex¬ 
periment  are  seen  to  agree  reasonably  well  with  the  exact  solution,  any 
apparent  discrepancy  being  possibly  due  to  the  intermittency  effects  men¬ 
tioned  earlier. 

Lastly,  the  distribution  of  the  ensemble  average  velocity  <U> 
at  prescribed  phase  angles  are  shown  in  figures  12  a,  b,  and  c,  for  a  few 
typical  points  in  the  oscillation  cycle,  for  each  of  the  runs.  The  distri¬ 
bution  of  the  time-mean  velocity,  0  is  also  shown,  in  each  case,  for  com¬ 
parison.  The  velocities  are  normalized  with  respect  to  the  time-mean 
velocity  at  the  center  line  (Umax)  corresponding  to  each  run.  The  distri¬ 
butions  at  the  higher  frequency  indicate  simple  displacement  but  very  little 
distortion  from  the  shape  of  the  time-mean  profile  in  the  core  region. 

Near  the  wall,  (in  the  Stokes  layer)  however,  the  profiles  cross  over  and 
suffer  some  distortion.  However,  no  reverse  flows  were  observed  in  the  pres¬ 
ent  experimrnts.  At  the  lower  frequency,  the  different  profiles  of  <U> 
show  distortion  with  respect  to  the  time  mean  profile  across  the  entire 
pipe  because  of  the  variation  of  both  the  amplitude  and  phase  angle  of  the 
periodic  velocity  continuously  across  the  pipe. 

C.  Effect  of  Flow  Modulation  on  Turbulence  and  Transition 

1.  General.  The  data  obtained  in  the  two  series  of  present 
experiments  provide  some  insight  into  the  effects  of  impressed  unsteadiness 
on  the  structure  of  turbulence  and  transition  in  pipe  flow.  For  convenience, 
the  two  series  of  experiments  will  be  discussed  one  by  one. 
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2.  Results  of  Series  1  experiments  (turbulent  unsteady  flow): 

The  distribution  of  the  rms  turbulent  intensity  u',  normalized  with  respect 
to  u*1  is  shown  in  figure  13.  Comparison  with  the  steady  flow  data  for 
9=0  shows  that  in  the  oscillatory  flow  there  is  a  slight  decrease  in  the 
maximum  turbulence  intensity  near  the  wall,  while  the  turbulence  intensity 
farther  away  from  the  wall  is  not  affected.  The  near-wall  region  where 
the  turbulent  intensity  is  affected  extends  approximately  over  25  per  cent 
of  the  pipe  radius  and  this  region  coincides  with  the  Stokes  layer  where 


significant  effects  of  flow  oscillation  are  observed  on  the  time-mean  and 
periodic  structures  of  the  flow  also.  This  result  is  in  agreement  with  the 
conclusions  of  Acharya  and  Reynolds  (1975)  that  the  structure  of  turbulence 
is  primarily  affected  within  the  Stokes  layer.  The  distributions  of  the 


ensemble  averaged  turbulence  intensity  Up(Q)  normalized  with  respect  to  the 
corresponding  shear  velocity  u*  (0)  are  also  shown  in  figure  13  for  a  few 
typical  fixed  phase  angles  of  the  oscillation  cycle.  Large  variations  in 


the  distributions  from  one  another  and  from  the  distribution  of  u'/u*  indi¬ 


cate  that  the  turbulence  structure  is  highly  disturbed  from  equilibrium. 


The  variations  in  the  relative  turbulent  intensity  (u'/u.  )  with  time  are 

P  P 

shown  in  figure  14  for  two  typical  locations  across  the  pipe.  These  figures 


suggest  that  the  flow  was  undergoing  rapid  structural  distortion  in  time. 
Actually,  the  wall  shear  stress  (u^  )  could  nearly  follow  the  changes  through 
the  discharge  cycle  while  the  turbulent  u-f luctuations  could  not.  This  is 
seen  clearly  from  figure  14  which  also  shows  the  variation  of  u^  and  the 

variation  of  u'  (not  normalized)  at  these  two  typical  points  in  the  pipe. 

P 

It  is  seen  that  u.  oscillates  with  a  slight  phase  lead  with  respect  to 
*P 

the  cross  sectional  average  periodic  velocity  Ump,  while  Up  remains 

^The  shear  velocity  u*  (0)  in  unsteady  flow  was  obtained  from  the  ensemble 
average  velocity  gradient  at  the  wall.  The  velocity  gradient  at  the  wall 
was  calculated  from  the  measured  <U(0)>  at  the  first  point  near  the  wall 
assuming  a  linear  velocity  distribution  in  the  viscous  sublayer  at  all 
instants.  This  procedure  was  found  to  be  generally  satisfactory  as  veri¬ 
fied  from  checks  made  in  steady  laminar  and  turbulent  flow.  In  these  cases, 
the  wall  shear  stress  obtained  by  this  method  compared  well  (within  3%) 
with  the  value  obtained  from  pressure  drop  measurements.  From  the  ensemble 
averaged  value  u*p,  the  time  averaged  value  u*  was  obtained  from  the  discrete 
analog  of  the  relation, 

i  2ir 

u?  =  d0 

*  2?r  J  *p 
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practically  constant  throughout  the  oscillation  cycle.  Such  a  freezing 

of  the  turbulent  normal  stress  can  be  expected  to  occur  when  the  flow  is 

subjected  to  rapid  strain  rates,  i.e.,  when  changes  occur  at  time  scales 

comparable  to  some  important  characteristic  time  scale  of  turbulence.  In 

a  wall-bounded  flow,  we  can  regard  the  turbulent  burst  frequency,  to 

determine  such  a  characteristic  time  scale  (at  least  near  the  wall).  In 

the  present  case  the  estimated  burst  frequency  is  about  2.3  Hz,  (using 
Um 

the  criterion  -r— r-  s  5  of  Rao  et  al  1971)  .  The  estimate  is  rough  but  is 
ttu 

sufficient  to  indicate  that  one  can  expect  to  find  interaction  between  the 
external  oscillation  and  the  internal  turbulence  structure  when  the  flow 
pulsates  at  a  frequency  of  1.75  Hz.  Stress-freezing  is  well  documented  in 
rapidly  accelerated  steady  flows  but  studies  of  rapidly  decelerated  steady 
flows  are  difficult  to  perform,  since  the  flow  in  such  cases  would  separate 
very  quickly.  One  of  the  interesting  features  of  the  present  experiments 
is  that  the  flow  experiences  large  pressure  gradients  of  alternating  sign 
but  yet  dees  not  actually  separate.  An  estimate  of  the  severity  of  the  pres¬ 
sure  gradient  in  the  present  case  can  be  made  by  calculating  the  value  of  the 
Clauser  pressure  gradient  parameter  3  defined  by  (neglecting  the  shear 
effects) , 


6  = 


D 


3P  .  D 


dU 


2pu?  9x  2u£  dt 


(23) 


Such  a  calculation  shows  that  6  varies  from  about  -40  to  about  +40  during 
a  cycle.  The  magnitudes  of  6  attained  during  the  cycle  are  comparable 
with  the  values  for  some  of  the  severe  adverse  pressure  gradient  steady 
flows  reported  in  the  Proceedings  of  the  Stanford  Conference  (Coles  and 
Hirst  1978).  Not  only  are  the  pressure  gradients  very  large  but  they  vary 
rapidly,  the  value  of  dg/dt  being  of  the  order  of  100  sec  1 ,  (based  on  a 
variation  in  8  from  40  to  -40  in  half  the  period  of  oscillation).  It  is 
thus  clear  that  the  turbulent  structure  will  be  in  a  highly  disturbed  state. 
Hence,  quasi-steady  turbulence  models  based  on  local  equilibrium  assump¬ 
tion  cannot  be  expected  to  describe  the  flow  at  such  oscillation  frequencies. 

The  quantity  which  is  of  greater  value  than  u'  in  understanding 
the  structure  of  turbulence  is  the  Reynolds  shear  stress  -  p  uv.  Unfortunately, 
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on  account  of  instrumentation  limitations,  it  was  not  possible  to  make 
direct  measurements  of  the  Reynolds  shear  stress.  However,  it  is  possible 
to  compute  it,  perhaps  with  some  loss  of  accuracy,  from  the  measured  wall 
shear  stress  and  velocity  distribution.  For  this  purpose,  we  write 
the  instantaneous  x-momentum  equation  for  pipe  flow  as 

3U  .  ..  3U  .  ..  3U  3P  ,  1  3 (t r) 

p  — •  +  pU  —  +  pV  —  =  -  —  + - e — *-  (24) 

v  3t  M  3x  M  3r  3x  r  3r 

where  V  is  the  instantaneous  velocity  in  the  radial  direction.  Performing 

an  ensemble  average  on  this  equation,  one  gets,  for  fully  developed  unsteady 
turbulent  flow, 

3<U> 


where  <t>  =  t, 


_  ^P>  1  3< (rr) > 

3t  3x  r  9r 

p<uv> 


(25) 

(26) 


laminar 

We  define  -  p<uv>  as  the  ensemble  averaged  Reynolds  shear  stress.  Mutli- 
plying  both  sides  of  the  equation  (25)  by  27rr  and  integrating  across  the 
pipe,  one  gets  (after  inserting  the  phase  average  notations  introduced 
earlier) 


3U  (6) 

m 


3t 


♦  4t»<9) 


3x 


D 


3P 


Eliminating  —  between  equations  (25)  and  (27)  one  gets 

aX 

r  2t  (6)r 

~(r , 6)  =  7  /  --  [<U(r,0)>  -  U  (0)]*r  dr  +  — 1 ^ - 

p  r  '  dt  m  D 


and 


-<uv> (r , 0)  =  ^  (r,6) 


3<U(r,9)? 


3r 


(27) 


(28) 


(29) 


The  ensemble  (phase)  averaged  Reynolds  shear  stress  -<uv>  was 
obtained  from  Equations  (28)  and  (29)  using  measured  velocity  and  wall 
shear  stress  data.  The  process  required  differentiation  with  respect 
to  time  and  integration  across  the  pipe  both  of  which  were  performed  numeri¬ 
cally  on  a  computer.  The  velocity-time  data  were  smoothed  in  order  to  re¬ 
duce  the  numerical  noise  in  differentiation.  The  procedure  was  found  to 
work  reasonably  well  because  the  velocity  values  were  available  at  very 
close  time  intervals  (50  values  per  cycle) ,  velocity  gradients  in  the 
radial  direction  were  not  severe  (due  to  the  low  Reynolds  number  of  the  flow) , 
sufficient  cycles  were  used  for  phase  averaging  (effectively  300)  and  lastly, 
the  three  terms  in  equation  (28)  were  all  of  similar  order  of  magnitude. 
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Variation  of  the  ensemble  averaged  Reynolds  shear  stress  during 
the  cycle  is  shown  in  figures  15  a,  b,  c  for  three  typical  locations  in  the 
pipe,  viz;  1,  2  and  10  mm  from  the  wall.  Although  there  is  some  scatter  in 
the  data,  the  results  are  still  good  enough  to  allow  one  to  arrive  at 
meaningful  conclusions  especially  in  the  region  of  flow  where  the  Rey¬ 
nolds  shear  stress  is  significant.  Very  near  the  wall,  <uv>  shows  a 
cylical  variation  but  lags  behind  (and  hence  behind  u^  also).  As 

the  distance  from  the  wall  increases,  the  amplitude  of  variation  of  <uv> 

decreases  until  it  becomes  almost  indistinguishable  from  scatter  beyond  the 
Stokes  layer  (as  seen  from  figure  15c).  The  well-defined  cylical  variation 
of  <uv>  in  the  inner  layer  (figures  15  a,  b)  is  in  contrast  to  the  frozen 

structure  of  u' .  This  suggests  that  the  modulation  in  <uv>  might  have 
P 

originated  from  modulations  in  the  ensemble  averaged  intensity,  v\  This 
is  possible  since  the  v-f luctuations  are  likely  to  be  of  finer  scale  than 
the  u-f luctuations  (by  a  factor  of  nearly  5  as  observed  by  Ramaprian,  1975) 
and  hence  can  respond  more  readily  to  the  impressed  oscillations.  Since 

the  time  scale  gets  larger  as  the  distance  from  the  wall  increases,  the 

response  of  <uv>  must  be  expected  to  diminish  eventually  in  the  core 
region.  The  nearly  frozen  structure  of  <uv>  in  the  core  region  (figure  15c) 
thus  indicates  rapid  strain  rates.  The  departure  from  local  equilibrium 
is  again  seen  from  the  variations  in  <uv>/u^  during  the  oscillation  cycle. 

Figure  16  a,  b,  c,  and  d  are  cross  plots  of  the  above  data 
and  show  some  typical  distributions  of  the  total,  laminar  and  ensemble 
averaged  Reynolds  shear  stresses,  across  the  pipe  for  fixed  phase  angles 
during  the  cycle.  The  distributions  are  normalized  using  the  relevant  inner 
layer  scales  of  velocity,  (uA  )  and  length,  (v/u.  ).  The  large  variations  in 

*  p 

the  distribution  of  the  stresses  from  one  point  in  the  cycle  to  the  other 
(seen  more  clearly  from  these  figures)  again  indicate  that  the  turbu¬ 
lence  structure  is  far  from  equilibrium  and  hence  cannot  be  described  by 


quasi-steady  models.  Figure  17  shows  the  time  averaged  distributions  of 
each  of  the  above  stresses  using  u^  as  the  scaling  velocity.  These  distri¬ 
butions  were  obtained  by  averaging  the  computed  distributions  at  50  phase 
positions  in  the  cycle.  It  is  seen  that  except  for  a  slight  discrepancy 
near  the  wall,  the  time-averaged  total  shear  stress  x  shows  the  expected 
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linear  distribution.  This  indicates  that  there  were  no  significant 
errors  in  the  numerical  procedure  used  for  calculating  the  shear  stresses 
from  the  experimental  data.  Some  of  the  important  features  pertaining  to 
figures  16  and  17  that  require  comments  are: 

(i)  The  laminar  shear  stress  is  significant  upto  a  consid¬ 
erable  distance  from  the  wall.  This  is  to  be  expected 
in  view  of  the  low  Reynolds  number  of  the  flow. 

(ii)  The  time  averaged  turbulent  shear  stress  reaches  a  maxi¬ 
mum  at  about  6  mm  (n=0.25)  from  the  wall  i.e.,  at  about 
the  same  location  where  u'  reaches  a  maximum.  It  is  impor¬ 
tant  to  note  that  this  is  within  the  Stokes  layer. 

(iii)  The  time  averaged  laminar  shear  stress  distribution  exhi¬ 
bits  a  local  minimum  near  the  wall.  This  is  due  to  the  pres¬ 
ence  of  a  point  of  inflection  in  the  time-mean  velocity  pro¬ 
file. 

(iv)  The  distribution  of  the  time-mean  Reynolds  shear  stress  also 
exhibits  a  kink  within  the  Stokes  layer.  Also,  it  is  ob¬ 
served  that,  in  the  region  very  close  to  the  wall,  the 
time-mean  Reynolds  shear  stress  attains  a  value  larger  than 
in  steady  flow. 

(v)  The  time-mean  distributions  of  the  stresses  are  not  signifi¬ 
cantly  different  from  the  corresponding  distributions  in 
steady  flow,  in  the  region  beyond  the  Stokes  layer.  However, 
by  affecting  the  distributions  within  the  Stokes  layer, 
the  imposed  unsteadiness  produces  an  increase  in  the  time- 
mean  velocity  gradient  at  the  wall  and  hence  increases  the 
wall  shear  stress.  This  is  seen  from  figures  14  a,  b, 
where  the  quasi-steady  time-mean  value,  u*  (obtained  from 
pressure  drop  measurements  at  the  mean  Reynolds  number  of 
2100)  is  shown  as  a  horizontal  dashed  line  against  the  dis¬ 


tribution  of  u.  in  unsteady  flow.  The  time-mean  value  of 
*P 

uA  is  clearly  larger  than  the  quasi-steady  mean  value  of  u^. 
Figure  18  shows  a  few  typical  distributions  of  the  ensemble 


average  velocity  <U>  in  the  inner  layer  co-ordinates  <0>/u.  vs  (D/2-r)u.  /v. 

*  P  *P 

The  time-mean  distribution  U/u^  vs  (D/2-r)u*/v  is  also  shown  in  the  figure. 


4 


28 


The  data  display  exactly  the  trend  that  one  would  expect  to  see.  In  the 

region  very  close  to  the  wall  y  <10,  the  ensemble  average  velocity  scales 

with  the  corresponding  u.  but  at  larger  distances  from  the  wall,  the  flow 

P 

structure  does  not  keep  in  step  with  uA  and  hence  the  data  at  different 
phase  angles  exhibit  different  distributions.  It  is  also  important  to  note 
that  there  is  a  net  difference  between  time  averaged  velocity  distribution 
in  the  unsteady  flow  (shown  by  open  circles)  and  the  velocity  distribution 
in  steady  flow  (shown  by  the  full  line).  This  difference  is  seen  more 
clearly  in  this  figure  than  in  figure  8,  discussed  earlier. 

The  quantity,  E^  =  j  {<x>  ~^-}2ir  r  dr  represent  the  rate  of  work 
done  per  unit  length  of  pipe  by  the  shear  stress  and  hence  its  average 
value  over  an  oscillation  cycle  represents  the  power  loss  in  the  pipe.  The 
term  normalized  with  respect  to  (pjP^irD/4)  is  shown  in  figure  19  as  a 
function  of  the  phase  angle.  The  corresponding  power  loss  in  a  hypothetical 
quasi-steady  turbulent  flow  can  be  obtained  for  the  same  Reynolds  number 
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modulation  using  the  Blasius  relation,  X  =  0.3165  Re  4.  Also,  for  the 
actual  transitional  quasi-steady  flow,  the  power  loss  for  the  same  Reynolds 
number  can  be  obtained  from  the  quasi-steady  pressure  drop  data  of  figure 
4  (or  alternately  the  friction  factor  data  of  figure  5).  Both  these  dis¬ 
tributions  are  also  shown  in  figure  19.  It  is  seen  that  flow  oscillation 
in  run  13  at  a  high  frequency  results  in  an  increase  in  average  power  loss 
when  compared  with  the  quasi-steady  transitional  flow.  This  is,  however, 
to  be  expected  since  the  unsteady  flow  is  fully  turbulent  during  the  entire 
cycle,  whereas  the  quasi-steady  flow  is  not.  However,  it  is  interesting 
to  observe  that  the  unsteady  flow  has  a  slightly  higher  average  power  loss 
when  compared  to  the  hypothetical  fully  turbulent  quasi-steady  flow  also. 


3.  Results  of  Series  2  experiments  (laminarized  unsteady  flow) 
During  the  second  series  of  experiments,  the  unsteady  flow  remained  nearly 
laminar.  Since  this  was  a  very  significant  observation,  it  was  carefully 
verified  to  insure  that  the  laminarization  was  not  due  to  any  obvious 
reasons  such  as  an  increase  in  fluid  viscosity,  reduction  in  flow  velocity, 
etc.  Starting  from  a  fully  turbulent  flow  at  the  fully  open  slot  position, 
the  flow  could  be  laminarized  by  just  rotating  the  sleeve  at  the  lowest  pos¬ 
sible  speed  (0.057  Hz).  The  photographs  in  figure  20  show  a  few  typical 
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oscilloscope  traces  obtained  from  the  LDA  at  several  oscillation  fre¬ 
quencies  and  include  data  from  two  radial  locations,  r  =  0  and  r  =  18. A 
mm.  In  each  photograph,  three  traces  are  shown:  (i)  fully  turbulent  steady 
flow  at  0=0  corresponding  to  a  Reynolds  number  of  about  2900,  (ii)  fully 
turbulent  laminar  flow  at  6=180°,  and  (iii)  fully  or  partially  laminarized 
oscillatory  flow.  The  photographs  generally  confirm  that  the  oscillatory 
flow  had  a  strong  tendency  towards  laininarization  under  the  operating  con¬ 
ditions  of  the  series  2  experiments. 

A  careful  study  of  the  traces  at  several  oscillation  frequencies 
indicated  that  laminarization  was  not  always  complete,  and  that  the  oscil¬ 
latory  flow  did,  in  fact,  become  turbulent  at  times  (as  seen  in  some  of  the 
photographs  in  figure  20).  The  intermittency  of  this  occurence  depended  on 
the  oscillation  frequency.  In  order  to  make  a  more  detailed  study  of  this 
phenomenon,  several  long-time  records  of  the  LDA  signal  were  obtanined 
using  a  strip  chart  recorder.  The  record  length  in  all  the  cases  was 
equal  to  300  cycles  or  5  minutes  whichever  was  longer.  Records  were  obtained 
for  r  =  0  and  r  =  18.4  mm.  Similar  records  were  obtained  in  quasi-steady 
flow  for  comparison.  A  few  typical  records  are  shown  in  figures  21  and 
22.  It  is  seen  that  the  flow  does  become  turbulent  several  times  during 
the  duration  of  the  record.  Estimating  the  intermittency  of  turbulence 
from  such  records  is  no  doubt  difficult  and  subjective;  particularly  when 
turbulent  and  oscillation  frequencies  are  not  well  separated.  The  following 
procedure  was  arbitrarily  selected  for  making  a  rough  estimate  of  the  inter¬ 
mittency.  The  procedure  is  based  on  the  observation  (and  assumption)  that 
whenever  the  flow  became  turbulent,  its  velocity  level  would  jump  upwards 
or  downward  from  the  laminar  value  depending  on  whether  the  point  under 
consideration  was  in  the  wall  region  or  on  the  centerline.  In  both  the 
cases,  a  distortion  would  be  observed  in  the  velocity  signal  recorded  on  the 
strip  chart.  The  dynamic  response  of  the  recorder  (5  Hz)  was  adequate 
for  indicating  this  distortion.  A  short  length  of  the  laminar  record 
traced  out  on  a  transparent  sheet  would  be  moved  over  the  record  and  portions 
of  the  record  that  did  not  coincide  with  the  laminar  record  would  be  marked 
out  as  turbulent  intervals.  The  intervals  marked  T  in  the  photographs  of 
figure  22  correspond  to  typical  turbulent  intervals.  The  intermittency 
factor,  y  was  defined  as  the  ratio  of  the  total  turbulent  interval  to  the 
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total  record  length.  In  the  case  of  steady  flow,  the  turbulent  intervals 
could  be  easily  recognized  without  much  ambiguity. 

The  intermittency  factor  is  plotted  in  figure  23  for  quaci-steady 
flow  as  a  function  of  quasi-steady  Reynolds  number.  It  is  seen  that  the 
mean  quasi-steady  flow  (6=90'’,  i.e.,  Re  s  2100)  has  an  intermittency  of 
about  0.2.  The  intermittency  factor  for  oscillatory  flow  is  shown  in 
figure  24  as  a  function  of  oscillation  frequency.  A  parameter  often  used 
in  characterizing  laminar  oscillatory  pipe  flows  is  .2  defined  as  [D/2/— ~ J . 

It  is  seen  that  fi  is  proportional  to  the  ratio  of  the  pipe  radius,  to  the 
thickness  of  the  Stokes  layer.  The  values  of  Q  are  also  shown  in  figure 
24.  The  data  for  unsteady  transitional  flow  show  some  interesting  trends. 

The  intermittency  is  very  low  at  either  end  of  the  frequency  range  studied 
(0=4-22).  However,  it  increases  fairly  significantly  and  appears  to  reach 
a  maximum  around  an  oscillation  frequency  of  about  0.4  Hz  (i2sl2)  .  This 
maximum  intermittency  factor  attained,  is  not  very  small,  and,  at  r=18.4  mm, 
is  indeed  comparable  to  that  of  the  mean  quasi-steady  flow.  Lastly,  the 
data  indicate  different  values  for  intermittency  near  the  wall  and  at  the 
center.  It  is  interesting  to  note  that  Yellin  (1966)  also  observed  that 
when  transition  occurred  in  pulsatile  flows,  the  disturbances  were  typically 
restricted  to  the  core.  However,  the  authors  have  no  explanation  at 
this  time  for  this  behavior  of  the  flow.  The  other  aspects  of  transition 
in  oscillatory  flows  are  examined  in  some  detail  below. 

Wygnanski  and  Champagne  (1973)  have  reported  detailed  study  of 
the  steady  transitional  pipe  flow.  They  divide  intermittent  transitional 
turbulence  into  two  categories  -  the  "puff"  type  and  the  "slug"  type.  Puffs 
are  caused  by  large  disturbances  at  the  inlet  region  of  the  pipe  and  are 
actually  remnants  of  a  partial  relaminarization  process.  Slugs,  on  the  other 
hand,  originate  from  the  instability  of  the  boundary  layer  and  represent 
various  stages  of  amplification  of  these  disturbances.  The  puffs  usuaLly 
occur  in  the  Reynolds  number  range  of  2000-2700  while  the  slugs  appear 
generally  at  very  much  larger  Reynolds  numbers.  The  intermittency  measure¬ 
ments  in  steady  flow  in  the  present  case  are  compared  in  figure  23  with  the 
puff  flow  data  from  Wyganski  and  Champagne.  The  qualitative  agreement  be¬ 
tween  the  two  sets  of  data  suggests  that  the  present  transitional  flow  can 
be  regarded  as  a  'puff  type  flow.  This  is  further  borne  out  by  the 


fact  that  the  turbulent  puffs  are  of  approximately  the  same  intensity  near 
the  wall  as  at  the  center  line  (as  seen  from  the  strip  chart  record),  unlike 
a  slug  type  flow  which  is  known  to  exhibit  an  increase  in  intensity  from 
the  centerline  to  the  wall.  It  may  be  noted  also  that  the  flow  is  fully 
turbulent  at  about  Re  =  2700. 

The  development  of  turbulence  through  a  puff-type  transition 
process  does  not  depend  on  shear  layer  instability  near  the  pipe  wall  since 
the  disturbance  is  provided  by  the  inlet  conditions  but  depends  on  the 
existence  of  a  Reynolds  number  above  a  critical  value  so  that  the  initial 
turbulence  can  be  sustained.  This  threshold  critical  value  is  about  2000 
and  transition  to  turbulence  can  not  occur  below  this  Reynolds  number. 

In  a  flow  subjected  to  a  favorable  pressure  gradient,  this  critical  Reynolds 
number  is  known  to  increase.  But,  in  an  adverse  pressure  gradient,  it  will 
not  decrease  significantly  below  the  threshold  value.  Hence,  if  the  Reynolds 
number  of  fluw  oscillates  around  2000,  the  net  effect  would  be  an  increase  in 
the  critical  Reynolds  number  and  hence  a  partial  or  complete  laminarization 
of  the  initial  turbulent  puffs.  The  increase  in  the  critical  Reynolds  number 
will  be  larger  at  larger  favorable  pressure  gradients  and  hence  one  would 
expect  to  find  a  greater  degree  of  laminarization  as  the  oscillation  fre¬ 
quency  increases.  This  is,  exactly  what  is  observed  in  figure  24  as  a  decrease 
in  the  intermittency  in  the  range  of  0.4  Hz  -  1.75  Hz  (ft=12-22)  with  the 
flow  being  almost  completely  laminar  at  1.75  Hz.  The  existence  of  a  maxi¬ 
mum  intermittency  (almost  equal  to  the  intermittency  in  the  mean  steady 
flow)  at  about  0.4  Hz  (fisl2)  appears  somewhat  to  contradict  the  above 
argument.  It  is,  however,  important  to  note  that  at  low  frequencies,  (0  - 
0.4  Hz),  (fi  =  0-12)  the  favorable  pressure  gradients  are  very  small  (at  the 
amplitude  of  modulation  employed)  and  consequently,  the  effect  on  the  criti¬ 
cal  Reynolds  number  is  negligibly  small  at  these  frequencies.  On  the  other 
hand,  the  extent  of  laminarization  of  the  puffs  depends  not  only  on  the 
critical  Reynolds  number  but  also  on  the  length  of  time  in  a  cycle  the  puffs 
are  exposed  to  a  Reynolds  number  lower  than  the  critical  Reynolds  number. 

At  low  frequencies  the  latter  effect  becomes  dominant  and  since  the  resi¬ 
dence  time  of  the  puffs  below  the  critical  Reynolds  number  increases  as 
the  oscillation  frequency  is  reduced,  the  Intermittency  of  puffs  would 
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decrease  with  frequency  in  this  range.  The  results  in  this  range  agree 
qualitatively  with  those  of  Sarpkaya  (1966)  whose  experiments  extended  over 
the  range  Q  =  4-7.8.  He  found  that  at  small  amplitudes  of  modulation 
(X  =  0.2  -  0.3)  there  was  only  a  slight  increase  in  critical  Reynolds  number 
above  the  steady  state  value  and  that  increase  in  the  oscillation  frequency 
in  this  range  reduced  the  rise  in  critical  Reynolds  number.  Sarpkaya  noticed 
negligible  effect  of  oscillation  on  the  transition  characteristics  at 
0  =  7.8  (which,  incidentally,  he  regarded  as  'rapid'  oscillation).  It  is 
interesting  to  note  that  at  about  0.4  Hz  corresponding  to  0^12  the  present 
unsteady  flow  exhibits  the  maximum  intermittency .  The  present  observations 
are  thus  in  qualitative  agreement  with  those  of  Sarpkaya.  There  is,  however, 
quantitative  disagreement  between  the  present  data  and  those  of  Sarpkaya 
with  respect  to  the  actual  magnitudes  of  the  critical  Reynolds  number  at 
the  various  values  of  SI.  This  is  due  largely  to  the  difference  between  the 
definitions  of  critical  Reynolds  number  used  in  the  two  cases.  In  the  present 
case,  the  critical  Reynolds  number  is  taken  to  be  the  Reynolds  number  at 
which  the  puffs  (external  disturbances)  disappear  whereas  Sarpkaya  defined 
it  as  the  Reynolds  number  at  which  the  external  disturbances  cease  to 
amplify.  It  is  thus  clear  that  the  present  definition  describes  the  lower 
bound  for  the  puff-type  transition  process. 

The  major  difficulty  with  studies  connected  with  puff-type  tran¬ 
sition  is  that  the  process  is  very  sensitive  to  the  nature  of  external 
disturbances  and  other  ambient  conditions.  It  is  thus  very  difficult  to 
find  repeatability  in  observations  over  a  period  of  time.  This  is  especially 
true  when  an  additional  factor,  namely  flow  oscillation,  is  introduced. 

This  has  been  clearly  demonstrated  by  the  two  series  of  experiments  reported 
in  this  paper.  The  two  series  gave  entirely  different  results  under  what 
apparently  appeared  to  be  identical  conditions.  It  is  worth  examining  the 
reasons  for  this  at  least  in  a  qualitative  way.  As  mentioned  already,  the 
mean  steady  flow  (0=90°)  was  fully  turbulent  in  the  Series  1  experiments. 

In  fact,  it  was  fully  turbulent  even  at  6=100°  (corresponding  to  Res2000) . 
Unfortunately,  intermittency  measurements  were  not  made  in  these  series. 
However,  one  can  visualize  (without  much  error)  the  intermittency  variation 
with  Reynolds  number  to  be  as  shown  in  figure  23.  It  is  very  clear  that 
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in  this  flow,  transition  occurred  within  a  very  narrow  range  of  Reynolds 
numbers,  unlike  in  the  second  series  of  experiments  in  which  it  was  spread 
over  a  much  wider  range  of  Reynolds  numbers.  The  reason  for  this  is  not 
known  with  certainty  but  a  reasonable  guess  is  as  follows.  The  first 
weries  of  experiments  were  conducted  in  summer  during  which  time,  the 
laboratory  temperature  was  about  30°C  while  the  fluid  temperature  remained 
at  about  26°C.  Since  the  pipe  was  not  insulated,  there  was  heat  transfer 
from  the  ambient  to  the  fluid.  This  would  aid  the  process  of  transition. 

On  the  other  hand,  the  second  series  of  experiments  were  conducted  in 
winter  when  the  laboratory  temperature  was  around  20°C  while  the  fluid 
temperature  still  remained  at  26°C.  There  was,  thus,  a  heat  transfer 
from  the  fluid  to  the  surroundings.  This  loss  of  energy  to  the  surroundings 
would  reduce  the  energy  available  for  producing  and  completing  the  transi¬ 
tion  and  hence  would  tend  to  stabilize  the  flow. 

The  two  intermittency  distributions  shown  in  figure  23  for  the 
two  series  can  be  used  to  explain  qualitatively  the  observed  difference  in 
response  of  the  flow  to  oscillate  at  a  high  frequency  in  the  two  cases. 

Flow  oscillation  at  a  high  frequency  would  raise  the  critical  Reynolds 
number  from  its  value  of  about  2000  at  quasi-steady  state  by  an  amount,  say, 
ARe.  This  is  equivalent  to  moving  the  operating  line,  i.e.,  mean  flow  Rey¬ 
nolds  number  of  the  unsteady  flow  to  the  left  by  ARe  in  figure  23. 

It  is  seen  that  in  the  series  2  experiments,  the  new  operating  line  cor¬ 
responds  to  a  larainarized  flow  (of  very  small  intermittency)  while  in  the 
series  1  experiments,  it  corresponds  to  fully  turbulent  flow  (of  intermittency 
1).  It  is  to  be  mentioned  here  that  laminarized  oscillatory  flow  can  be  ob¬ 
served  only  under  very  restricted  conditions.  These  include  puff-type  transi¬ 
tion  (i.e.,  transition  brought  about  by  inlet  disturbances),  a  relatively 
small  intermittency  of  puffs  at  the  mean  Reynolds  number  and  either  strong 
pressure  gradient  fluctuations  (high  frequency,  large  amplitude  oscillation) 
or  large  puff-residence  time  (long  enough  pipe).  These  conditions  are  often 
found  in  the  pulsatile  flow  of  blood  in  the  mammalian  aorta  and  hence  lamin- 
arization  of  the  flow  can  be  expected  to  occur  in  such  cases.  This  is,  in 
fact,  corroborated  by  some  of  the  recent  in-vivo  aorta  experiments  mentioned 
at  the  beginning  of  this  paper. 


IV.  CONCLUSIONS 


The  present  study  has  lead  to  the  following  conclusions: 

(i)  Periodic  oscillation  of  discharge  tends  to  increase  the 
critical  Reynolds  number  of  puff-type  transitional  pipe  flow.  Under  cer¬ 
tain  conditions,  the  transitional  flow  may  be  laminarized  on  periodic 
oscillation.  For  a  given  amplitude  of  flow  modulation  the  extent  of 
laminarization  depends  on  factors  such  as  the  intermittency  of  puffs  in  the 
quasi-steady  mean  flow,  the  oscillation  frequency  and  the  residence  time 

of  the  puffs  in  the  pipe. 

(ii)  The  laminarized  periodic  flow  behaves  very  much  like 
laminar  periodic  flow.  For  example,  the  time-mean  flow  properties  re¬ 
main  unchanged  from  those  of  quasi-steady  mean  flow  and  the  phase  lag 
and  amplitude  of  the  periodic  velocity  component  depend  strongly  on  the 
Strouhal  number. 

(iii)  When  the  oscillatory  flow  is  fully  turbulent,  its  time 
mean  and  periodic  structure  qualitatively  resemble  those  of  oscillatory 
laminar  flow  at  the  same  Strouhal  number.  However,  the  behavior  of  the 
oscillatory  turbulent  flow  is  also  influenced  by  an  additional  parameter, 
namely  the  ratio  of  the  oscillation  frequency  to  some  characteristic 
frequency  of  turbulence.  When  this  ratio  is  of  the  order  unity,  the  oscil¬ 
lations  interact  with  the  turbulent  structure.  Important  differences 

can  be  observed  between  laminar  and  turbulent  flows  at  such  oscillation 
frequencies.  For  example,  the  time-mean  velocity  profile  in  the  oscillatory 
flow  exhibits  a  point  of  inflection  near  the  wall,  and  the  time  mean  wall 
shear  stress  and  power  loss  increase  from  their  quasi-steady  values.  Also, 
the  periodic  velocity  component  exhibits  an  overshoot  in  the  Stokes  layer, 
the  magnitude  of  the  overshoot  being  larger  than  in  laminar  oscillatory 
flows  at  the  same  Strouhal  number. 

(iv)  At  the  interactive  frequency  of  scillation,  mentioned  above, 
the  ensemble  averaged  turbulence  intensity  is  frozen  everywhere  in  the  pipe. 
The  ensemble  averaged  Reynolds  shear  stress  is  able  to  follow  the  oscillation 
cycle  (with  some  lag)  only  very  close  to  the  wall.  However,  beyond  the 
Stokes  layer,  it  is  also  frozen  at  some  average  value.  The  stress  freezing 
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is  brought  about  by  the  large  and  rapidly  varying  strain  rates.  The 
ensemble  averaged  Reynolds  stresses  as  well  as  the  ensemble  averaged 
velocity  do  not  scale  with  the  corresponding  ensemble  averaged  wall 
shear  stress  indicating  significant  departures  from  local  structural 
equilibrium. 
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Figure  1.  Schematic  layout  of  the  experimental  apparatus 
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Oscilloscope  traces  of  velocity  signals 
From  I.DA  (r  =0)  .  The  smooth  horizontal  trace 
corresponds  to  steady  laminar  flow  (G=180°) ; 
the  hashed  horizontal  trace  corresponds  to 
steady  turbulent  flow  (0=0°) ;  the  smooth  sin 
usoidal  trace  corresponds  to  laminarized 
oscillatory  flow 


Figure  7.  Velocity  distribution  in  steady  flow  in  wall-layer  coordinates 
9=0  (Maximum  slot  opening) 


Figure  8.  Distribution  of  the  time-mean  velocity  across  the  pipe 
in  unsteady  flow. 

A,  steady  laminar  flow  (9  =180°) ;  -6-,  steady  turbulent  flow  (9=0°), 
“•7  unsteady  turbulent  flow  at  f=1.75  Hz  (Run  13); 

®,  unsteady  laminarized  flow  at  f=1.75  Hz  (Run  23), 
k,  unsteady  laminarized  flow  at  f=0.057  Hz  (Run  24) ; 

- ,  theoretical  parabolic  profile  for  laminar  flow 


0.0  90.0  183.0  270.0  3S0.0  9$0.0  S90.0 

9 (degrees) 

Figure  9 (b) .  Variation  of  the  periodic  component 
angle  (Run  23) 


Figure  10(c).  Distribution  of  the  periodic  component  of  the  velocity  across  the  pipe  at  fixed  phase  angle 
in  oscillating  laminar  flow;  f =0.057  Hz  (Hun  24) 


Fiqure  11.  Distribution  of  the  amplitude  of  the  periodic  component  of  velocity  across  the  pipe. 

f,  Run  13;  ®,  Run  23;  A  ,  Run  2A; - ,  Laminar  theory  from  Uchida  (1956) 


Figure  12 (a) .  Distribution  of  the  ensemble  average  velocity  across  the 
pipe  in  oscillatory  turbulent  flow,  f=1.75  Hz  (Run  13) 


symbol,  phase  angle  ©(degrees):  •,  45°;  C,  135°;  8,  180°; 

®,  225°;  0,  270°;  A,  315°,  A,  360°  time-mean  velocity 


Figure  12 (b) .  Distribution  of  the  ensemble  average  velocity  across  the 
pipe  in  oscillatory  laminarized  flow,  f=1.75  Hz  (Run  23) 
symbol,  phase  angle  ©(degrees):  0,  45°;  C,  135°;  9,  180° 
0,  225  ;  #,  315°;  A,  360";  — -  time-mean  velocity 


Qqqqooooqgqo 


Figure  15(a).  Variation  of  the  Reynolds  shear  stress  with  ohase  angle  in  oscillatory 
flow,  distance:  1  mm  from  the  wall;  f=1.75  Hz  (Run  13) 


9 (degrees) 

Time  variation  of  the  Reynolds  shear  stress  with  phase  angle  in 
oscillatory  flow,  distance:  10  mm  from  the  wall;  f=1.75  Hz  (Run  13) 
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Distribution  of  the  time-averaged  total,  laminar  and 
Reynolds  shear  stresses  across  the  pipe . 

Steady  turbulent  flow  at  a  Reynolds  number  of  2870  (6=0°) 
I  ,  total  shear  stress;  Q  ,  Reynolds  shear  stress; 

0,  laminar  shear  stress 

Unsteady  turbulent  flow,  f=1.75  Hz  (Run  13) 

•  ,  total  shear  stress;  A,  Reynolds  shear  stress; 

A,  laminar  shear  stress 


Figure  20.  Oscilloscope  traces  of  velocity  signals  from  LDA 
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Figure  20.  Oscilloscope  traces  of  velocity  signals  from  LDA  from 
experiment  series  2.  r=18.4  mm 
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Figure  21.  Typical  strip  chart  records  of  the  LDA  output  signal  in  steady  flow 
from  experimental  series  2 
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APPENDIX  A 


Description  of  the  LDA  System  at  IIHR  and  sample 
calculations  of  velocity  from  LDA  signals 


I 
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Description  of  the  IIHR  LDA  System 


A.  The  LDA  System 

The  simple  channel  LDA  system  at  IIHR  is  shown  schematically 
in  figure  A. 1.  The  system  consists  of  a  SPECTRAPHYSICS  5  mw  He-Ne 
laser,  ill'minating  TSI  (Thermo  System  Inc.)  optics.  The  optics  consists 
of  a  beam  splitter  which  splits  the  laser  beam  into  two  beams  each  of 
half  the  original  strength  and  spaced  50  mm  apart.  Each  of  the  beams  is 
then  passed  through  an  accousto-optic  Bragg  cell  which  shifts  the  fre¬ 
quency  of  the  laser  light  by  a  preset  amount.  The  Bragg-cells  and  the 
associated  electronics  were  supplied  by  Opto-Elektronische-Instrumente. 

It's  possible  to  shift  the  frequency  of  each  beam  independently  by  7 
selectable  discrete  frequencies  ranging  from  38.0  MHz  to  44.0  MHz. 

By  shifting  the  two  beams  by  different  frequencies  one  can  obtain  a  rela¬ 
tive  frequency  shift  between  the  two  outcoming  beams.  In  the  present 
case,  a  relative  frequency  shift  of  0.3  MHz  was  used.  The  dual  beams  from 
the  Bragg  cells  are  focussed  by  the  focussing  lens  of  104.1  mm  focal 
length  at  the  point  (in  reality,  a  finite  sampling  volume)  in  the  fluid 
where  the  velocity  is  to  be  measured.  The  light  scattered  by  a  particle 
which  passes  through  this  volume  is  picked  up  by  the  receiving  optics 
shown  in  the  figure.  The  receiving  optics  consists  of  a  pair  of  lenses  and 
a  prism  that  focusses  the  scattered  light  on  to  a  photodiode.  The  photodiode 
converts  the  optical  signal  to  an  electrical  signal  which  is  then  amplified. 

In  this  particular  arrangement,  known  as  the  dual  beam  system, 
the  scattered  light  is  essentially  a  mixture  of  the  light  scattered  by 
the  particle  from  the  two  incident  beams.  Let  be  the  original  laser 
beam  frequency,  and  be  the  frequency  shifts  caused  by  the  Bragg 
cell  on  the  two  incident  beams,  n^  and  n^»  the  unit  vectors  in  the 
directions  of  the  incident  beams,  r  the  unit  vector  in  the  direction  of 
collection  (see  figure  A-l).  Then,  for  an  idealized  case  of  scattering 
occurring  due  to  a  single  particle  travelling  with  a  velocity  V  as  shown, 
assuming  uniform  plane  waves,  one  can  write 


»i  •  +  VB1 + 


(A.  1) 
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r-n 


v2  =  Vo  +  VB2  +  tV"T 


i2 1 


(A. 2) 


where  and  \>2  are  the  frequencies  of  scattered  light  from  the  two 
beams,  and  the  bracketed  terms  represent  the  Doppler  shift  frequency  in 
each  case.  The  light  picked  up  by  the  photodiode  is  a  mixture  of  the 
above  two.  From  the  principle  of  hetrodyning,  it  can  be  shown  that 
the  resulting  light  will  contain  the  mixture  of  the  sum  and  difference 
of  the  above  two  frequencies.  Of  these,  the  difference  equal  to 


(v 


B1 


'  VB2) 


+  V*  [ 


n . „-n  , 
i2  il. 


(A. 3) 


is  the  frequency  of  interest.  This  frequency  is  detected  in  the  output 
signal  by  a  frequency  tracker.  A  TSI  model  1090  tracker  is  used  for 
this  purpose.  The  tracker  output,  which  is  a  linear  function  of  the  in¬ 
stantaneous  velocity  of  the  flow  is  low-pass  filtered  at  150  Hz  to  remove 
unwanted  noise  and  then  connected  to  the  IIHR  Data  Acquisition  System.  The 
low-pass  filtering  frequency  of  150  Hz  was  found  to  be  high  enough  to  include 
all  significant  turbulent  frequencies  in  the  flow.  It  is  seen  from  Eq. 

A. 3  that  the  output  signal  frequency  is  independent  of  the  collection 
direction  r  (through  the  intensity  of  the  collected  light  does  depend  on 
r) .  It  is  also  seen  that  at  zero  velocity,  the  output  frequency  will  be 
v  ^  -  Vg2  i.e.,  equal  to  the  imposed  relative  frequency  shift. 

It  is  appropriate  to  mention  here  briefly  some  of  the  problems 
encountered  in  LDA  measurements.  One  of  them,  viz.,  directional  ambiguity 
has  been  overcome  in  the  present  studies  by  using  frequency  shifting  tech¬ 
nique.  Secondly,  particle  size  and  distribution  introduce  complexities, 
not  encountered  in  the  simple  single  particle  analysis.  The  particle  size 
influences  scattered  light  distribution  and  intensity,  the  ability  of  a  par¬ 
ticle  to  follow  the  flow  and  the  signal  to  noise  ratio.  Also,  the  higher 
the  particle  concentration,  the  more  will  be  the  background  noise  due  to 
phase  differences  between  particles  arriving  at  the  sampling  volume. 

Low  particle  concentrations,  on  the  other  hand,  limit  the  velocity  versus 
time  information  as  well  as  generally  increase  the  complexity  of  data  processing 
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due  to  discontinuity  in  the  data,  or  "signal  drop  out".  Signal  drop 
out  can  also  be  caused  by  large  variations  in  velocity  with  time  as  can 
occur  in  highly  turbulent  flows,  very  low  velocities  and  in  oscillatory 
flows.  In  these  cases,  the  accompanying  frequency  changes  are  too  large 
to  be  tracked  by  the  tracker.  The  use  of  frequency  shifting,  in  fact,  has 
the  added  advantage  of  bringing  the  range  of  frequency  variations  into  the 
most  efficient  region  of  tracker  operation.  Ideally,  the  particles  should 
have  a  high  refractive  index,  be  uniform  in  size,  and  be  as  large  as  pos¬ 
sible  (up  to  5y)  while  still  following  the  flow.  The  concentration  should 
be  low  enough  to  minimize  Doppler  ambiguity  while  high  enough  to  obtain  ( 
adequate  time  information  on  flow  characteristics.  To  approach  this  ideal 
seeding  is  generally  required.  However,  naturally  occurring  particulates 
in  the  fluid  may  often  be  adequate  to  give  satisfactory  results,  as  in  the 
present  study.  Thirdly,  the  spatial  resolution  of  LDA  is  limited  by  the 
requirement  of  a  finite  sampling  volume.  With  reference  to  figure  A-l  the 
dimensions  and  D^  of  the  sampling  volume  depend  on  the  intersection 
angle  <$>  of  the  beams  in  the  fluid  and  the  diameter  of  the  beam.  In  the  pre¬ 
sent  experiments  the  dimensions  Lm  and  are  1.1  mm  and  0.7  mm  respectively. 
The  dimension  is  very  much  larger  than  the  diameter  of  a  hot-wire  (2.5  pm) 
but  comparable  to  the  diameter  of  a  total  head  tube.  It  is,  however,  to  be 

noted  that  since  oil  was  used  in  the  present  studies,  the  dimension  L  is 

m 

still  less  than  the  thickness  (*1.5  mm)  of  the  viscous  sublayer.  The 
calculations  required  for  obtaining  the  value  of  velocity  from  the  LDA 
output  are  described  in  Appendix  B. 

B.  The  Traverse  Mechanism  for  LDA 

A  simple  one  dimensional  traverse  mechanism  was  built  for 
traversing  the  LDA  optics  along  a  horizontal  diameter  of  the  test  section. 

The  construction  of  this  traverse  can  be  seen  from  figure  A-2.  The 
upper  alluminium  channel  slides  over  the  bottom  alluminium  channel.  The 
latter  rests  on  a  mild  steel  framework  and  can  be  adjusted  to  be  level 
by  using  the  four  levelling  screws.  The  upper  channel  is  moved  longi¬ 
tudinally  by  a  spindle  and  has  a  travel  of  50  mm.  A  dial  indicator  of  0.01 
mm  least  count  is  mounted  in  such  a  way  as  to  measure  the  displacement  of 
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the  upper  channel.  The  laser,  the  entire  optics  and  the  photodiode  are 
all  mounted  on  wooden  bases  fixed  to  the  upper  alluminium  channel  and  thus 
move  together  as  one  unit. 

Because  of  the  refraction  effect,  through  the  tube  wall  and 
the  liquid,  the  actual  displacement  of  the  focal  point  (sampling  volume) 
is  not  equal  to  the  physical  displacement  of  the  traverse  (apparent  dis¬ 
placement).  Knowing  the  thickness  of  the  pipe  wall  and  the  refractive 
indices  of  the  pipe  wall  and  the  oil,  one  can  calculate  the  actual  displace¬ 
ment  of  the  focal  point  for  a  given  apparent  displacement.  It  was,  how¬ 
ever,  felt  that  a  practical  calibration  would  be  more  desirable,  as  well  as 
concenient,  in  the  present  case.  A  temporary  arrangement  consisting  of 
an  L-shaped  point  gage  traversed  by  a  precision  micrometer  (0.025  mm 
least  count)  was  used  for  this  purpose.  The  point  gage  was  traversed 
through  a  hole  in  the  tube  wall,  located  at  such  a  position  that  the  tip  of 
the  gage  travelled  along  the  LDA  axis.  The  location  of  the  focal  point 
could  be  easily  determined  from  observing  the  light  scattered  by  the  tip  of 
the  gage  as  it  moved  along  the  axis  of  the  LDA.  The  focal  point  was  given 
a  series  of  prescribed  apparent  displacements  as  measured  by  the  dial  indi¬ 
cator  and  the  corresponding  actual  displacements  as  measured  by  the  micro¬ 
meter  were  obtained.  The  resulting  calibration  curve  is  shown  in  A. 3  and 
is  seen  to  be  linear  (as  expected) .  This  result  was  used  in  all  subsequent 
experiments  to  obtain  the  actual  displacement  of  the  focal  point.  Theoreti¬ 
cally,  the  uncertainty  involved  in  the  calibration  was  equal  to  the  dimen¬ 
sion  "Lm"  of  the  sampling  volume.  In  practice,  however,  it  was  much  less 
than  this  since  it  was  possible  to  locate  the  two  "extremes"  of  this  region 
and  set  the  point  gage  at  the  middle  of  these  two  positions.  The  finite 
length  of  the  sampling  volume  did  present  some  difficulty  in  fixing  the 
"zero"  distance  from  the  wall  during  the  velocity  traverses.  This  diffi¬ 
culty  was  overcome  by  matching  the  velocity  profile  in  steady  laminar  flow 
with  the  theoretical  parabolic  distribution.  The  "zero"  so  determined 
was  used  for  the  turbulent  flow  traverse  also.  In  the  unsteady  flow, 
the  location  of  the  "zero"  position  was  taken  to  be  the  point  at  which 
the  oscillation  amplitude  (as  observed  on  the  oscilloscope)  would  suddenly 
go  to  zero.  This  was  a  very  well  defined  point  and  could  be  located  with 
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excellent  repeatability. 


C.  Sample  Calculation  of  Velocity  from  the  LDA  Signal 


If  fp  is  the  Doppler  frequency  shift,  from  Eq.  (A.  3) 


fD  XL*  (nil~ni2) 


From  Fig.  A.l 

v.  c  nV 

hence,  tQ  =  f- 


Ini2"nill 


ni2~nil 


2  sin  4> / 2 

2  ^  sin  <+>/ 2 

AL 


—r  -r- 

where  U  is  the  velocity  component  normal  to  the  bisector  of  n_^  and  n^* 

The  value  of  n  is  taken  as  unity  in  air  and  the  angle  as  tan--*-  (d/2L,) 
where  d  is  the  original  beam  spacing  introduced  by  the  beam  splitter  and 
L  is  the  focal  length  of  the  focussing  lens.  It  can  be  easily  shown 
that  even  when  the  beams  pass  through  media  of  different  refractive  indices 
(in  the  present  experiment,  for  example,  they  pass  through  the  plexiglass 
tube  wall  and  then  through  the  oil),  the  produce  (n  •  sin  <f)/2)  can  still 
be  calculated  using  the  above  values  for  n  and  $.  Hence  Eq.  (A. 3)  can  be 


written  as 


2U  sin  [Tan  (d/2L) 

X. 


U  =  fD^2  sin  [tan-1 (d/2L) K  fD 

For  the  present  LDA  system,  X^=  632.8  nm,  L  =  104.1  mm  and  d  =  50  mm. 


Hence, 


K  =  1.313  m/s/MHz 


If  the  relative  frequency  shift  introduced  by  the  Bragg  cells  is  f  and 

the  measured  frequency  is  f 

U  =  (f  -f  )  K 
m  s 


As  already  mentioned,  the  value  of  f  in  the  present  experiments  was  0.3 
Hz.  The  tracker  output  voltage  E  is  proportional  to  f  or  E  =  cf^. 


Hence, 


v  v  (£-En ) K 

0  =  (-  -  fs)K  (E-Cfs)  -  = 
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where  Eq  is  the  tracker  output  voltage  at  zero  velocity.  In  he  present 
experiments  c  was  1  volt/MHz  1  0.4%.  Thus, 

U  -  (E  -  0.3)  x  1.313  m/s 


DESCRIPTION  OF  THE  1IHR  LDA  IN  DUAL  BEAM,  FORWARD  SCATTER  MODE 


DETAILS  OF  THE  LDA  TRAVERSE 


apparent  displacement 
actual  displacement 


FIG. A3.  CALIBRATION  OF  THE  LDA  TRAVERSE 


83 


APPENDIX  B 

Experimental  Data 


AD-A081  180  IOWA  INST  OF  HYDRAULIC  RESEARCH  IOWA  CITY  F/6  20/4 

EXPERIMENTS  ON  TRANSITIONAL  OSCILLATORY  PIPE  FLOW*(U) 

AU6  79  B  R  RAMAPRI AN»  S  W  TU  DAAG-29-79-S-0017 


UNCLASSIFIED  IIHR-221  NL 


END 


OATI 

--/•LMED 

3-80 


DOC 


84 


54  2.033 


60 


66  |  2.575 


72 


78  3.184 


84 


90  3-715 


96 


102  4.304 


108 


114  4.854 


120 


132 


1381  5.965 


144 


150  6.486 


156 


162  7.075 


168  7.390 


Uave  :( 


TABLE  1. 

STEADY  PRESSURE  DROP 


-5  2 

v=  1.384  x  10  m/s 
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_2 
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3.854 
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2.552 


3.018 


4.098 
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3.339 
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5.399 


6.726 


7.236 


|  7.864 
m/s  )  , 


6.777 


6.547  6.282 


7.726 
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7.153  6.649 


6.943 


_ 7.795  |  7.189 

dp/dx  :  (nevrtan/m2/m) 


7.789 


6.326 


8.821 


3.567 


3.915 


4.262 


4.608 


4.956 


5.253 


5.649 


5.945 


6.076 
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5.958 

5.840 
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7.831 

8.424 


9.564 


MEASUREMENTS  IN  STEADY  FLOW 
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TABLE  3. 

MEASUREMENTS  IN  UNSTEADY  FLOW  (  RUN  23  ) 


°inax=1.086in/s,  (Uqn)pmax=  0.43  m/s  *  u«  =0-°3^  m/s 
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APPENDIX  C 


Assembler  Language  Program  for  Sampling  and 
Processing  Periodic  Turbulent  Flows  on  the  IIHR 
IBM/1800  System 
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//  JOB  0000011111 
//  ASM  COMPT 


•LIST 

• 

PROGRAM. 

COMPT 

• 

THIS  PROGRAM  COMPUTES  THE  PERIODIC  IIEAI.' 

• 

AND  THE 

RMS  OF  THE 

RANDOM  COMPONENT 

• 

OF  THE 

SIGNALS  FOR 

3  CHANNELS  FROM  THE  DATA 

• 

STORED 

ON  DISK  BY 

THE  PROCCGS  PROGRAM  AGS V 3 . 

* 

PROGRAM 

MED  BY  A. MUELLER.  APRIL  1977 

0000 

01 

0C000152 

.  START  XlO 

L  LI  TON 

0002 

01 

C40006F9 

LP 

L  3TDT1 *1 

SAVE  SECTOR  ADDRESS. 

0004 

01 

04  0000  7A 

STO 

L  SECSV 

0006 

01 

C4000370 

LD 

L  DE320 

SET  UP  WORD  COUNT  TO  OCTAI 

N 

0008 

01 

04000103 

STO 

L  IIIFOM 

EXPERIMENTAL  PARAMETERS. 

OOOA 

20 

04262495 

LICF 

DISKN 

READ  F  IL2  0  FROM  DISK 

0006 

0 

1000 

DC. 

yiooo 

WHICH  CONTAINS  THE 

OOOC 

1 

0163 

DC 

INFOM 

NUMBER  OF  POINTS  PER  WAVE, 

0000 

0 

0000 

DC 

/oooo 

THE  NUMBER  OF  WAVES  SAMPL- 

ooof: 

20 

04262495 

LIBF 

DISKN 

ED  AND  THE  RUN  NUMBER. 

OOOF 

0 

0100 

DC 

/0100 

0010 

1 

0163 

DC 

INFOM 

0011 

0 

73FC 

MDX 

•  -4 

0012 

01 

C400016S 

LOO  PI  LO 

L  WAVES 

SET  UP  COUNTER  FOP  NUf'CFR 

0011. 

01 

D40O015A 

STO 

L  f.'UMCR 

OF  AVERAGES  COMPUTED. 

0016 

01 

C40001C  7 

LD 

L  POUTS 

COMPUTE  TWICE  THE  NUMBER 

0018 

Q 

1001 

SLA 

1 _ 

OF  POINTS. 

0019 

01 

D4000158 

STO 

L  TNPTS 

001B 

0 

10A0 

SLT 

32 

CLFAR  STORAGE  FOP. 

001C 

01 

65800158 

LDX 

11  TNPTS 

ACCUMULATION  OF  SUMS. 

001E 

01 

DD0002A6 

CLEAR  STD 

1.1  AREAC-2 

AREAC  :  USED  FOR  ACCUMULA- 

0020 

0 

71FE  . 

ilDX 

1  -2 

TION  OF  SUMS. 

0021 

Q 

70FC 

MDX 

CLEAR 

0022 

00 

6500FDC0 

LDX 

LI  -576 

CLEAR  AREAD. 

0024 

01 

D50005B  5 

CLER2  STO 

Ll  AREAD*  5  76 

0026 

0 

7101 

MDX 

1  1 

0027 

0 

70FC 

MDX 

CLER2 

002B 

01 

94000158 

S 

L  TNPTS 

OBTAIN  MINUS  TWICE  TifE 

002A 

01 

D4000158 

STO 

L  TNPTS 

NUMBER  OF  POINTS.  TNPTS 

002C 

01 

C400007C 

LD 

L  XR4 

CONTAINS  -2 •PONTS. 

_QQ2E 

01 

C400007CL 

CMP 

L  TWO  . 

CHECK  FOR  LOOP.  SET  UP  OF 

0030 

0 

7002 

MDX 

USKIP 

ADDRESSES. 

0031  0 

7010 

MDX 

SKIP 

0032 

0 

70 1A 

MDX 

SKIP 

0033 

0 

CO  2  F 

NS KIP  LD 

AODDl ♦ 1 

0034 

01 

94000158 

s 

L  TNPTS 

MODIFY  ADDRESSES  FOR 

0036 

0 

D02C  __ 

STO 

ADDDl ♦  1 

ACCUMULATION  AMR 

0037 

0 

D02H 

STO 

AODDl  +  3 

MAf.'l  PULA  TION  OF  DATA. 

_QQ38 

01 

D 4 000099 

STO 

L  LDDD1 ♦ 1 

003A 

01 

C400011C 

LD 

L  ADDD2^1 

SET  UP  ACDRESSES . 

00 

01 

940001 58 

s 

L  TNPTS 

DO  FOR  F 1 RST  LOOP  ONLY. 

003E 

01 

D400011C 

STO 

L  ADDD2»1 

0040 

01 

D400011E 

STO  L  ADDO  2*  3 

0042 

0 

1010 

SLA 

16 

ZERO  A  REG. 

0043 

01 

94000167 

__  s 

L  POUTS 

COMPUTE  MINUS  THE 

004S 

01 

D4000159 

STO 

L  MMPTS 

NUMBER  OF  POINTS. 

0047 

01 

C4 000167 

LD 

L  PONTS 

0049 

01 

84000114 

A 

L  SURCP+1 

004B 

01 

04000114 

STO 

. L  SUBCP*1 

0040 

01 

C40001 57 

SKIP  LO 

L  WDCMT 

SET  UP  WORD  COUNT  FOR 

004F 

01 

D4000GF8 

STO 

L  STOT1 

DISK  READ. 

0051 

20 

04202495 

REAUD  LIDF 

DISKN 

READ  IN  THE  FIRST  SET  OF 

0052 

0 

1000 

DC 

71000 

DATA  STORED  ON  DISK  III 

< 


I 


0053 

o8i% 

_ 005& 

0057 

0058 

1 

So 

0 

1 

0 

U6F  8 

DC 

nr 

ST0T1 
/  0000 

04202495 

0100 

LI  BF 
DC 

DISKM 

/0100 

0GF8 

70  FC 

DC 

MDX 

STDT1 

•  -4 

0059 

01 

65800158 

AGAIN  LDX 

"FT 

tnpTs 

load  Index  reBTstfRS  T6 

0050 

01 

GC800159 

LOX. 

1  2 

MNPTS 

BEGIN  SIGNAL  AVERAGING. 

0050 

01 

C4000GF A 

LOAD1  LD 

i 

5TDT1 ♦ 2 

_ D05F 

0 

D019 

$TO 

SAVEl 

0060 

0 

C01  8 

ID 

SAVFl 

0061 

0 

189  0 

SRT 

16 

COMPUTE  SIGflAL  AVERAGED 

0062 

01 

8D0002A8 

A0DD1  AD 

Ll 

AREAC 

RESULTS. 

.  0064 

QL 

DD0QQ2A8 

STp 

Ll 

ARE  AC 

0066 

01 

7401005E 

MDX 

L 

LOAD  1 *1/1 

MODIFY  ADORFSS. 

0068 

0 

C0F5 

LD_ 

LOAC1 ♦ 1 

CHECK  FOR  Eflp  OF  DATA  BLOCK. 

0069 

01 

94 0001 50 

s 

L 

ADCOl 

0Q6B 

-ILL 

04000157 

CMP 

L 

l.'DCNT 

006D 

0 

4  OOF 

BSI 

MORDT 

MORE  OATA  REQUIRED. 

0066 

0 

7001 

MDX 

MCOMP 

COMPUTATIONS  CONTINUE. 

006F 

0 

4  00D 

BSI 

MORDT 

MORE  DATA  REQUIRED. 

0070 

0 

7201 

MCOMP  MDX 

2 

1 

DECREMENT  XR2 . 

0071 

0 

1000 

NOP 

0072 

0 

7102 

MDX 

1 

♦2 

DECREMENT  XR1. 

0073 

0 

70E9 

MDX 

LOAD1 

GO  BACK  FOR  .'.ORE  DATA. 

0074 

_01 

74FF015A 

MDX 

_L  . 

NUMCR, -1 

DECREMENT  NUMBER  OF  l.AVES. 

00  76 

0 

7QE2 

MDX 

AGAIN 

GET  ANOTHER  WAVE  ^ET . 

0077 

0 

701C 

MDX 

COMPT 

GO  TO  COMPUTATION  SCHEME. 

0078 

0 

0000 

TEMP  DC 

•“A* 

TEMPORARY  STORAGE. 

0079 

_Q _ 

0000 

SAVEl  DC 

*  -  * 

TEMPORARY  STORAGE. 

007A 

0 

0000 

SECS V  DC 

•  -  * 

SAVE  SECTOR  ADDRESS  FOR 

007S 

Q 

QQQ2 

Tv;a  dc 

2 

RAVI  DATA. 

007C 

0 

0003 

XR4  DC 

3 

INDEX  REGISTER  TOR  LOOP. 

*  DISK  DATA  SUBROUTINE 

TO  OBTAIN  MORE  DATA.  * 

0070 

0 

0000 

MORDT  DC 

•  -  • 

SUBROUTINE  TC  GET 

007E 

01 

C4000158 

LD 

L 

ADCOl 

MORE  DATA  F ROM  DISK. 

0080 

01 

D400005E 

STO 

L 

LOADl^l 

RESET  ADDRESSES  OE  LOADING 

0082  01 

D40001 1 2 

STO 

L 

LOAD  2 ♦ 1 

OPERATION. 

0084 

01 

C40006F9 

LD 

L 

STDTl ♦ 1 

updatf  SECTOR  ADDRESS. 

0086 

_Q1_ 

8  4  00  01 6  2 

A 

L 

DECO 

0088 

01 

D40006F9 

STO 

L 

STDTl ♦ 1 

_ PQ8A 

20 

04262495 

- _ LLBf 

DIStti  _ 

CALL  TO  CISK.F0R 

*080 

0 

1000 

DC 

/iooo 

MORE  DATA. 

_ QQ8C  1 

06F  S 

_  J)C  _ 

STDTl 

008D 

0 

0000 

DC 

/oooo 

OGEE 

20 

Q42G2495 

LIGF 

disk:; 

008F 

0 

0100 

DC 

/0100 

Q090 

1 

06F8 

DC 

STDTl 

0091 

0 

70FC 

MDX 

•  -4 

_0092 

01 

4C80007D 

_  .  BSC 

j 

MORDT 

END  SUBROUTINE. 

0094 

01 

GSRomsa 

COr.PT  LDX 

11 

TNPTS 

LOAD  INDEX  REGISTERS 

0096 

01 

66800159 

LDX 

1  2 

MNPTS 

TO  SCALE  THE  OATA. 

0098 

01 

CD0002A3 

lddui  lpd 

Ll 

AREA  C 

009A 

01 

AC00016S 

D 

L 

WAVES 

009C 

01 

04000508 

STOR2  STO 

L 

CB'^CD*  2 

UNSCALED  SIGNAL  AVERAE,  0 

0096 

01 

7401009D 

MDX 

L 

STOR2M,  1 

RESULT  STORED  IN  F  1  L2 

QQAO 

Jl_ 

7102 

MDX 

♦2 

RESULT  STORED  BACK  INTO 

00A1 

0 

iooo 

NOP 

AREAC. 

00A2 

0 

7201 

MDX 

2 

♦1 

00A3 

0 

70F4 

MDX 

LDDOl 

00A4 

01 

74FF007C 

MDX 

L 

XR4, -1 

DECREMENT  OF  INDEX  REG. 

00A6 

0 

7001 

MDX 

LUCKV 

FOR  LOOP  COUNTER. 

00A7 

0 

7010 

MDX 

LOOP4 

ill  in  in  mmiitummi 
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00A.8 

01 

C400007C 

lucky  in 

L 

XR4 

BBac  H1  '?8i,ru,u 

/$£ 

Ihnpit 

OOAD 

0 

70on 

MDX 

LOOP  3 

OOAF 

0 

7000 

MOX 

L00P2 

_  QQAF 

Q1 

C<4000F.7C 

I00P2  LD 

L 

C  TD  T  2 ♦ 1 

SET  UP  SECTOP  ADDRESS  FOR 

OOB  1 

01 

U4000CF9 

STO 

L 

S  TDT1 ♦ 1 

CHANNEL  2. 

OOB  3 

01 

C4  00011*0 

LD 

L 

ADCOl 

CHANNEL 2  A  RESET  LOAD 

OOB  8 

01 

D400005E 

STO 

L 

L0AD1 ♦ 1 

INSTRUCTION. 

0007 

01 

4C000012 

BSC 

L 

L0UP1 

CO  BACK  TO  OPFRATE  Oil 

0009 

01 

C4000E7E 

L00P3  LD 

L 

STDT3M 

CHANNEL  2  S  GET  SUM.  OF  CHAN. 3. 

OOBC 

01 

DL0C0GF9 

STO 

|f 

S  T  0  T 1  ♦  1 

OOBD 

01 

C4000153 

LD 

L 

ADCOl 

QQBF 

01 

D4000Q5E 

STO 

L 

L0AD1 ♦ 1 

00C1 

01 

4CG00012 

nsc 

L 

L00P1 

GO  BACK  TO  OPERATE  ON 

Q0C3 

01 

C4U00Q7A 

LOOP 4  LD 

L 

SECSV 

CHANNEL  3. 

00C5 

01 

D40006F9 

STO 

L 

STDT1*1 

SAVE  SECTOR  ADDRESS  FOR 

QQC7 

Q1 

74030070 

MDX 

_L_ 

XR4.3 

CHANNEL1. 

00C9 

0 

CO  05 

LO 

CONST 

SAVE  SECTOR  ADDRESS  FOR 

.  OOCA 

01 

DL0005C6 

_ STO 

L 

CBJCD 

CBJCD.  _ 

oocc 

20 

04  2G2495 

LI  OF 

D  1  S  KN 

STORE  UNSCALED  SIGNAL 

QQCD 

0 

3000 

.  ..  DC  . .. 

/  3000 

AVERAGED  DATA  (AREAE) 

OOCE 

1 

05B6 

oc 

CB  .CD 

INTO  DISK. 

.  QQCF 

Q 

0140 

rnriST  nc 

320 

OODO 

20 

042G2495 

lief 

D  1  S  KN 

0001 

0 

0100 

..  .  DC  . 

/0100 

00D2 

1 

0  53  6 

DC 

CBJCD 

0003 

c 

70  FC 

MDX 

•  -4 

0004 

01 

C4000168 

INDEX  LP 

L 

WAVES 

0006 

01 

0  4  Q  Q  0 1 5  A _ 

STO 

_L_ 

CUMPR 

0008 

01 

C4C0007C 

LD 

L 

XR4 

LOAD  INDEX  REGISTER  TO 

OOPA 

01 

C 4  00  0  0  73 

CMP  L 

TWO 

CHECK  THE  LCCP  STEP. 

oooc 

0 

70  28 

MDX 

READF 

0000 

0 

7014 

MOX 

LOOPB 

CODE 

01 

C4000E73 

!  00 PA  LD 

L 

STDT2 ♦ 1 

LOAD  A  STORE  THE  DATA  OF 

00  E  0 

01 

D4Q00GF9 

STO 

L 

STDT1 ♦ 1 

CHANNFL  2  G. CHANGE  ADDRESS 

00E2 

01 

C4000114 

LO 

L 

SUBCP^i 

oo  r.  4 

01 

840001G7 

A 

L 

PONTS 

OOFG 

01 

04000114 

STO 

L 

SUBCP+1 

OOFS 

01 

C400015B 

LD 

L 

ADCOl 

00EA 

01 

04000112 

STO 

L 

LOAD 2*1 

00EC 

0 

CO  2  F 

LO 

A0DD2 ♦ 1 

MODIFY  ADDRESS  FOR 

OOEO 

01 

94000158 

S 

L 

TNPTS 

ACCUMULATION  -  MANIPULATION 

QOEF 

0 

DO  2C 

STO 

A  D  D  D  2  ♦  1 

OF  DATA. 

OOF  0 

0 

DO  2  D 

STO 

ADDD2»3 

00F1 

0 

7013 

MDX 

READF 

OOF  2 

01 

C4000E7E 

100P0  LD 

L 

S  TO  T3  ♦  1 

LOAD  A  STORF  THE  DATA  OF 

OOF  4 

01 

P40006F9 

STO 

L 

STDT1 *1 

CHANNEL  3  fi  CHANGE  ADDRESS. 

00F6 

01 

C4000114 

LD 

L 

SUBCP*1 

OOF  8 

01 

C40001G7 

A 

L 

PONTS 

OOFA 

01 

04000114 

STO 

L 

SUCCP*1 

OOFC 

01 

C,  4000150 

LD 

L 

AOCOl 

OOFE 

01 

D4000112 

STO 

L 

LOAD  2 ♦ 1 

0100 

0 

C01 3 

LD 

ADDD2*1 

MODIFY  ADDRESS  FOR 

0101 

01 

94000158 

S 

L 

TNPTS 

ACCUMULATION  5  MAN  1  PULAT  1  01.' 

0103 

0 

DCl  8 

STO 

APDP  ? ♦ 1 

OF  DATA. 

0104 

0 

DO  1  3 

STO 

A  D  D  0  2  ♦  3 

OIOS 

20 

042G2495 

READF  LIRF 

DISKN 

READ  OUT  THE  SET  OF  DATA 

010G 

0 

1000 

DC 

/  1000 

STORED  ON  DISK. 

0107 

1 

06F  8 

DC 

S  TOT  1 

0108 

0 

0000 

DC 

/OOOO 

0109 

20 

042G2495 

LI  DF 

DISKII 

010A 

0 

0100 

DC 

/0100 

010D 

1 

0GF8 

DC 

STDT1 

010C 

0 

70  FC 

MOX 

•  -4 

010D 

01 

G58001S8 

AGAN1  inx 

1  1 

TNPTS 

LOAD  INDEX  REGISTER  TO 

100 


010F 

—0111 

0113 

0115 

0117 

0119 

01 

01 

01 

01 

01 

01 

6G800159 
C4D00CF A 
9GGU05i<3 
A4U0015C 
D400U078 
A4000U78 

LOAD  2 
SUCCP 

LDX 

LD 

r 

•J 

M 

12 

-L 

L2 

L 

MNPTS 
STDT1 ♦ 2 
AREAE-1 
SCALE 

BEGIN  SIGNAL  AVERAGING. 

LOAD  FROM  DISK  AREA.  _ 

SUBTRACT  UNSCALED  MEAN 

RAW  DATA  u  MEAN  VALUE. 

STO 

M 

L 

L 

TEMP 

TEMP 

oTTfT 

-ST 

“SOTTO  376' 

A5U02 

Ad 

TT 

Ahead -l 

0110 

01 

DD000374 

STD 

LI 

AREAD-1 

01 1 F 

01 

74010112 

MDX 

L 

LOAD 2 ♦l t 1 

MODIFY  ADDRESS. 

0121 

0 

CO  FO 

LD 

LOAD2*l 

CHECK  FOR  END  OF  DATA  BLOCK 

0122 

01 

9400015  P* 

S 

L 

ADCOl 

012  6 

01 

P,  4000157 

CMP 

I. 

V»DCHT 
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PROGRAM  :  CllVUT 
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Y)  THROUGH  THE  SCALING  PROCEDURE. 
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1 

0346 

DC 

CB3CD 

0016 

0 

70FC 

MDX 

*-4 

0017 

01 

C4000208 

LD 

L 

PORTS 

0019 

0 

1001 

SLA 

1 

001 A 

01 

D 4 0001 40 

STO 

L 

TNPTS 

001C 

0 

1010 

SLA 

16 

SET  UP  ZERO  IN  A-REG. 

001D 

01 

94000208 

S 

L 

POUTS 

SUBTRACT  NO.  OF  PTS . 

001F 

01 

D  4  0  0  0 1 3  F 

STO 

L 

MNPTS 

MINUS  NO.  OF  PONTS 

0021 

01 

65S0013F 

LDX 

11 

MNPTS 

DO  LOOP  FOR  CHANNEL1 

0023 

01 

C 4 000348 

LOOP1 

LD 

L 

AREAE-1 

LOAD  (MEAN  ♦  PERI) 

0025 

01 

84000142 

A 

L 

ZEP.Ol 

ADD  UNSCALED  ZERO  VOLY 

0027 

01 

A4  0001 43 

M 

L 

SCALE 

SCALING 

0029 

01 

D 4 000144 

ARRAG 

STO 

L 

AREAG 

00  2D 

01 

74010024 

MDX 

L 

LOCPl+1,1 

0020 

01 

74  0100  2A 

MDX 

L 

ARRAG+1, 1 

002F 

0 

7101 

MDX 

1 

1 

INCREASE  INDEX  REG. 1  BY-! 

0030 

0 

70F2 

MDX 

LOOP1 

GO  TO  LOOP1 

0031 

0 

C0F2 

LD 

LOOP1+1 

0032 

01 

94000208 

S 

L 

POUTS 

0034 

0 

DOEF 

STO 

LOOP1+1 

0035 

0 

C0F4 

LD 

ARRAG+1 

0030 

01 

94000208 

S 

L 

PONTS 

>1038 

0 

D0F1 

STO 

ARRAG+1 

0039 

01 

C 4 000488 

LD 

L 

A  ORE 

CALCUL.  ADRESS  3.CH. 

003G 

01 

84000140 

A 

L 

TNPTS 

00  3D 

01 

9400013A 

S 

L 

DEC1 

003F 

0 

D003 

STO 

LOOP2+1 

0040 

01 

6680013F 

LDX 

1  2 

MNPTS 

DO  LOOP  FOR  CHANNEL  3 

0042 

00 

C4  00000  0 

LOOP2 

LD 

1. 

LOAD  MEAN  ♦  PERIODIC 

0044 

01 

A 4 000143 

M 

L 

SCALE 

0046 

01 

D40001A4 

ARRAU 

STD 

L 

AREAH 

0048 

01 

74010043 

MDX 

1. 

LOOP2+1, 1 

004A 

01 

74010047 

MDX 

L 

ARRAH+1, 1 

004C 

0 

7201 

MDX 

2 

1 

RAISE  INDEX  REG. 2  BY  1 

004D 

0 

70F  4 

MDX 

LOOP2 

108 


004E 

0 

C0F  7 

LD 

ARRAU 

oour 

01 

94000208 

S 

L 

POUTS 

0051 

0 

DCF  5 

ST0 

ARRAII*1 

0052 

20 

23A17171 

L 1  OF 

TYPE1 

PRINT  RUNNO 

0055 

1 

020A 

DC 

CLICK 

00514 

0 

0001 

DC 

1 

0055 

01 

XlD 

T~ 

■HWWJM 

0057 

01 

E40000E7 

AIJD 

L 

MASC2 

TEST  FOR  SUITCII  2 

0050 

0 

4820 

BSC 

Z 

SKIP  IF  DES  2  IS  NOT  SET 

005A 

0 

7000 

MDX 

SKIPP 

SKIPP  OUTPUT  CHANNEL  1 

0053 

20 

26665A15 

LI6F 

URTYN 

WRITE  MESSAGE  FOR  SEDIMENT 

035C 

0 

2001 

DC 

/  2  00 1 

0050 

1 

0489 

DC 

UMES1-1 

005E 

0 

0000 

DC 

/oooo 

005F 

20 

26663A15 

LI  OF 

URTYN 

0060 

0 

0001 

DC 

/0001 

0061 

0 

70FD 

MDX 

«-  3 

0002 

01 

04000208 

LD 

L 

PONTS 

00G4 

0 

D00  2 

STO 

MOP 

0065 

20 

25A17171 

ur>F 

TYPE1 

TYPE  OUT  SED.  CONCENTRATION 

0066 

1 

014  4 

DC 

AREAG 

IN  PPM  UNIT 

0067 

0 

0000 

NOP 

DC 

#  -  * 

0063 

20 

26663A15 

SKIPP 

LI  BF 

URTYN 

WRITE  MESSAGE  FOR  V-COMP 

0069 

0 

2001 

DC 

/  2  00 1 

0C6A 

1 

049A 

DC 

DUES  2-1 

006D 

0 

0000 

DC 

/oooo 

0060 

20 

26663A1 5 

LIPF 

URTYN 

006D 

0 

0001 

DC 

/  0001 

0Q6E 

0 

70FD 

MDX 

«-3 

_ QOGF 

01 

C4000208 

LD 

L 

PONTS 

0071 

0 

D002 

STO 

NNOP 

0072 

20 

23A17171 

LI  DF 

TY  PEI 

TYPE  OUT  V-COMP. 

0073 

1 

01A4 

DC 

AREAH 

IN  100Q*(FT/S)  UNIT 

00  74 

0 

0000 

NNOP 

DC 

*-A* 

0075 

01 

0CC000E4 

XIO 

L 

IOCCA 

READ  DATA  ENTRY  SUITCII 

0077 

01 

F.40000F.6 

ADD 

L 

MAS  C 

0079 

0 

4806 

CSC 

♦ 

007A 

0 

7067 

MDX 

OUT 

0070 

20 

03059115 

LIDF 

CARDN 

FEEDS  1  CARD 

007C 

0 

3000  , 

DC 

/  3000 

0070 

0 

0000 

DC 

0 

007F 

20 

03059115 

LI  OF 

CARON 

00  7F 

0 

0000 

DC 

/OOOO 

OOSO 

0 

70FD 

MDX 

.  *-?  _ 

0081 

01 

C400020A 

LD 

L 

CLICK 

PUNCH  RUNNO 

0083 

20 

02255103 

El  BF 

Bitmc 

003  4 

1 

00E9 

DC 

OUTPT 

WfT-T'l 

01 

LD 

L 

l.'AVES 

PUNCH  NO  OF  WAVES 

0087 

20 

02255103 

L 1  BF 

BINDC 

0083 

J. 

006F 

DC 

OUT  PT  +  6 

00  3  9 

01 

C400013D 

LD 

L 

DEC12 

008D 

0 

D050 

STO 

AREA 

008C 

20 

03059115 

L 1  BF 

CARDN 

MiMM 

0 

FTiTTtH 

DC  .. 

/  2000 

0086 

1 

00ES 

DC 

AREA 

008F 

0 

0000 

DC 

0 

0090 

20 

03059115 

LI  BF 

CARDN 

0091 

0 

0000 

DC 

/OOOO 

0092 

0 

70FD 

MDX 

*-  3 

■rrm 

01 

FFTTTrai 

L 

DEC  72 

PUNCH  MEAN  AND 

0095 

0 

D052 

STO 

AREA 

PERIODIC  ON  CARDS 

009G 

31 

65300208 

LDX 

11 

PONTS 

0098 

0 

6200 

pimcn 

LDX 

2 

12 

12  VALUES  PER  CARD 

0099 

01 

C40001A4 

CONV 

LD 

l 

AREAH 

0090 

20 

02255103 

LIUF 

BINDC 

0090 

1 

OOE9 

OTP 

DC 

OUTPT 

109 


009D 

01 

740C009C 

MDX 

L 

OTP,  6 

FIELD  FOR  NEXT  VALUE 

j.'.pX 

l 

rnuv* l f 

1  NEXT  VALUE . 

iiDX 

2 

00A2 

0 

70  FG 

MDX 

CON  V 

OOA3 

20 

03059115 

L1BF 

CARDN 

_ noAii 

0 

2000 

PC 

/  2  000 

00A5 

1 

00E8 

PC 

AREA 

0  0  A  0 

0 

0000 

DC 

0 

00 A  7 

20 

03059115 

LIOF 

CARDN 

QQAS 

0 

0000 

oc 

/  0000 

00A9 

0 

70FD 

IIDX 

*-  3 

0  0AA 

01  74080090 

MDX 

L 

OTP, -72 

RESTORES  1 .  OUTPUT  AD. 

OOAC 

0 

71F4 

MDX 

1 

-12 

OOAD 

0 

70EA 

MDX 

PUNCH 

OOAE 

01 

C4000208 

LD 

L 

pouts 

PRINTOUT  OF  REST 

OODQ 

0 

1890 

SRT 

16 

oon 

01 

AC00013D 

D 

L 

DEC12 

000  3 

0 

1090 

SIT 

1G 

0004 

01 

04000139 

CMP 

L 

ZERO 

0000 

0 

7002 

MUX 

GOON 

0007 

0 

7024 

MDX 

RSTOR 

0008 

0 

7023 

MDX 

RSTOR 

0009 

01 

114  00  0 141 

gooii 

STO 

L 

REST 

0000 

01 

A4000130 

M 

L 

DEC6 

OOCD, 

,  0 

1090 

SLT 

16 

OOOE 

01 

D40000E8 

STO 

L 

AREA 

OOCO 

0 

C0D9 

LD 

CONV+l 

SET  ADRESS  IN  AREAH 

OOCI 

0 

DOOG 

STO 

CONVV+1 

OOG2 

01 

84000141 

A 

L 

REST 

_ on  rii 

0 

nnns 

STO 

cor;v+i 

0CC5 

01 

GG800141 

LDX 

1  2 

REST 

00C7 

00 

C4000000 

cotiw 

LD 

L 

*  -  * 

0  0C9 

20 

02255103 

LI  BF 

BINOC 

OOCA 

1 

00E9 

OTTP 

BC 

OUTPT 

OOCB 

01 

740GOOCA 

MDX 

L 

CTTP, 6 

nor  n 

01 

IiDX 

L 

CONVV+1 

,  1 

0  0  OF 

0 

72  r  F 

MDX 

2 

-1 

0000 

0 

70  F  6 

MDX 

COM  VV 

00D1 

20 

03059115 

LIOF 

CARDN 

0002 

0 

2000 

DC 

/  2000 

0003 

1 

O0E8 

DC 

AREA 

0004 

DC 

0 

00D5 

20 

03059115 

LIOF 

CARDN 

OODG 

0 

0000 

DC  . 

/  0000 

0007 

0 

70F0 

MDX 

»-  3 

0008 

0 

C0F1 

LD 

CTTP 

RESTORES  OTTP 

00D9 

01 

940000E8 

S 

L 

AREA 

0 

doff 

STO 

OTTP 

OODC 

01 

C400009A 

RSTOR 

ld" 

L 

CONV+l 

RESTORES  ADR.  IN  CONV+l 

oonc . 

,  01 

94000208 

s 

L 

POUTS 

OOF.O 

01 

U40UQ09A 

STO 

L 

CONV+l 

00E2 

30 

O59C9SC0 

..  OUT 

CALL 

EXIT 

00E4 

0000 

OSS 

E 

0 

0  OF  4 

0 

IOCCA 

DC 

7  0000 

READ  DATA  ENTRY  SW. 

OOF.  5 

0 

0  740 

DC 

/  0  74  0 

00E6 

0 

0001 

I  IAS  C 

DC  . 

/0001 

OOF  7 

0 

0002 

MA3C2 

DC 

/0002 

MASK  FOR  DES  2 

OGES 

0 

0000 

AREA 

DC 

*-A  * 

OUTPUT  AREA  FOR  CARDS 

0QF9 

0050 

OUT  PT 

RSS 

so 

01  39 

0 

noon _ 

ZERO 

_ac _ 

0 

013A 

0 

0001 

0EC1 

DC 

1 

01313 

Q 

000G 

DECG 

DC  „ 

6 

013C 

0 

noGS 

DEC10 

DC 

104 

01 3D 

0 

oooc 

0EC12 

DC 

12 

013E 

0 

0048 

DEC72 

DC 

72 

013F 

0 

0000 

mnpts 

DC 

#-# 

0140 
0141 
0142 
0143 
0144 
0 1 A  4 

0 

0 

0 

0 

0000 

Qfion 

TNPTS  DC 

REST  DC  - 

*- * 

0000 

2717 

ZEROl  DC 

SCALE  DC- 

0  ZERO  SHIFT  CHAN. 1 

10007 

0060 

0060 

AREAG  BSS 
AREAII  C6S 

96  PPM  FOR  SED.  COIJC. 

96  1 0 0 0 » < FT/S )FOR  V-COMP. 

C  2  04 

31 

0C253CC0 

INFOM  DSA 

F  1  L20 

0207 

0 

0000 

DELAY  DC 

*  -  * 

0  20S 

0 

0000 

POUTS  DC 

*  -  * 

0209 

0 

0000 

WAVES  DC 

*  -  * 

0  2  0A 

0 

0000 

CLICK  DC 

#-#  . 

0  2  or. 

Htmni 

96 

0260 

0006 

AREAF  DSS 

219 

0346 

31 

00253C67 

C3jlCD  DSA 

FIL27 

0349 

013F 

AREAE  BSS 

319 

043S 

1 

0349 

ADRE  DC 

AREAE 

0439 

0 

0010 

DC 

DfiESA-D.MESl 

0  4  8  A 

0009 

DIIESl  DUES 

'  2  RPRES  SURE  ' 

0  481- 

dmes 

(MEAN  AND  PERIODIC)  'E 

0499 

0 

DC 

/  8105 

049A 

VtYnjrH 

DUESA  BES 

0 

049A 

0 

wrrtm 

DC 

DMESD-DMES  2 

0490 

■tYi  vm 

dues 2  d;;es 

1 2RVELOC 1 TY  COMPONENT* 

■lit  n 

dues 

(MEAN  AND  PERIODIC) 'E 

04AF 

0 

8105 

DC 

/  8105 

0400 

0000. 

DM  ESP.  RES 

0 

0  400 

0000 

END 

START 

1 

'!0_ . 

ERRORS.  10 

ABOVE  ASSEMBLY. 

cr:  vrt 
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//  job  oooooillll 

//  ASM  CONVT _ 

•LIST 


*  PROGRAM  :  CONVT  • 

*  THIS  PROGRAM  CONVERTS  VOLTAGE  UNIT  OF  RAHPOM _ *_ 

*  FLUCTUATIONS  F  OP,  PRINTOUT.  AND  PUNCHES  CARDS.  * 

»  DES 1  UP  :  CARDS  ARE  PUNCHED _  _ *_± 

*  DES2  UP  :  PRINTOUT  CHANNEL  1  IS  SUPRESSED  ~  * 

«  PROGRAi'lMEO  BY  A.  MUELLER _ MARCH  19  77 _ *_ 

*************************************************** 


0000 

01 

C400022C 

START  LD 

L 

DEC9 

LOAD  WORD  COUNT  ON  FIL20 

0002 

01 

1)4000234 

STO 

L 

1 11 F  Of  1 

0004 

20 

04262495 

LIBF 

DISKN 

GET  INFORMATION 

0005 

0 

1000 

DC 

/1000 

0006 

1 

0234 

DC 

INFOM 

0007 

0 

0000 

DC 

/OOOO 

OOOS 

20 

04262495 

LI  CF 

D  1  SKN 

0009 

0 

0100 

DC 

/0100 

OOOA 

1 

0234 

DC 

INFOM 

ooor. 

0 

70FC 

MDX 

•  -4 

oooc 

01 

C  4  0  0  0  4  8  7 

LD 

L 

WORDC 

LOAD  WORD  COUNT  ON  FIL22 

ooor 

01 

D 4 000242 

STO 

L 

CPVST 

0010 

20 

042G2495 

LI  CF 

DISKN 

GET  SQUARE  VALUE  OF 

0011 

0 

1000 

DC 

/  1000 

RANDOM  FLUCTUATION 

0012 

1 

0242 

DC 

CPVST 

0013 

0 

0000 

DC 

/OOOO 

0C1  4 

20 

04262495 

LIBF 

DISKN 

0015 

0 

0100 

DC 

/0100 

0016 

1 

0242 

DC 

CPVST 

0017 

0 

70FC 

MDX 

»-4 

0012 

0 

1010 

SLA 

16 

SET  A  REGISTER  ZERO 

001  J 

01 

94000238 

S 

L 

PONTS 

ooir. 

01 

D400022D 

STD 

L 

MNPTS 

MINUS  NO.  OF  PONTS 

001D 

0 

1001 

SLA 

1 

: i ULTI PLY  BY  TWS 

ooir 

01 

D400022E 

STO 

L 

TNPTS 

-2*  NO.  OF  PONTS 

0020 

0 

1001 

SLA 

1 

0021 

01 

D400022F 

STO 

L 

FNPTS 

-4*  NO.  OF  PONTS 

0023 

01 

C4000239 

LD 

L 

WAVES 

0025 

0 

1090 

SRT 

16 

SET  UP  DOUBLE  PRECISION 

00  26 

30 

0510G4E3 

FLOAT  NO.  OF  WAVES 

0028 

20 

058A3580 

LI  BF 

ESTO 

STORE  FLOATED  WAVES 

0029 

1 

0230 

DC 

DVSER 

002A 

20 

23A17171 

LIBF 

TY  PEI 

PRINT  RUNNO 

002B 

1 

0  2  3A 

DC 

CLICK 

002C 

0 

0001 

DC 

1 

0020 

01 

0CG00112 

XIO 

L 

IOCCA 

READ  DATA  ENTRY  SWITCH 

002F 

01 

E 4 000115 

AND 

L 

MASC2 

TEST  FOR  SWITCH  2 

0031 

0 

4820 

BSC 

Z 

SKIP  IF  DES  2  IS  NOT  SET 

0032 

0 

703  7 

MDX 

SKIPP 

SKIPP  OUTPUT  CHANNEL  1 

0033 

20 

2G663A15 

LIBF 

WRTYN 

CARRAGE  RETURN 

0034 

0 

2001 

DC 

/  2 001 

0035 

1 

0488 

DC 

0036 

•0 

0000 

DC 

0 

0037 

20 

2C6C3A15 

LIBF 

WRTYN 

0038 

0 

0001 

DC 

/0001 

0039 

0 

70FD 

MDX 

*-S 

20 

LIBF 

WRTYN 

0036 

0 

2001 

DC 

/  2001 

00  3C 

1 

0494 

DC 

DUES  1-1 

HEADING  FOR  CHANNEL  1 

00  3D 

0 

0000 

DC 

0 

003E 

20 

26663A15 

LIBF 

WRTYN 

003F 

0 

0001 

DC 

/0001 

0040 

1 

003E 

DC 

•  -3 

0041 

88H1 

00  44 
0045 
0  01*6 

20 

? 

0 

20 

0 

2GGG3A15 

UUF 

WKTYN 
/  7001 

CARRAGE  RETURN 

0000 

Be 

DC 

CAR  1 G 

0 

20G63A15 

0001 

LI  OF 
DC 

WRTYN 

/0001 

• 

00  4  7 

0 

70  FD 

MDk 

*-  3 

OO'tfi 

01 

6580022E 

I.OX 

11 

TfJPTS 

SET  UP  INDEX  REG.  1 

0  0 1»  A 

01 

CC000244 

L00P1 

LDD 

L 

AREAD-1 

LOOP  FOR  CHANNEL  1 

0  01*C 

30 

051064E3 

CALL 

EDFLT 

FLOAT  THE  VALUE 

00  l*E 

20 

05100940 

LIBF 

ED  1 V 

DEV  IDE  BY  NO.  OF  WAVES 

0  0 1*  F 

1 

0230 

DC 

DVSER 

0050 

30 

05898640 

CALL 

ESQR 

TAKE  SQUARE  ROOT 

0052 

20 

091S99CO 

LI  OF 

!  F  1  X 

CONVERT  TO  INTEGER 

0053 

01 

D400010C 

ADO 

STO 

L 

AREAG 

0055 

01 

74010054 

MDX 

L 

AD0*1, 1 

0057 

01 

740700 49 

IIDX 

L 

L00P1 ♦ 1 , 2 

0050 

0 

7102 

L0OPC 

nnx 

1 

2 

00  5  A 

0 

70F.F 

MDX 

L00P1 

005B 

01 

C4000238 

LD 

1. 

PONTS 

0  0  50 

0 

0002 

STO 

UPON 

00  5  E 

20 

2  3A 1  7 1  7 1 

LIBF 

TYPE1 

PRINTOUT  1. CHANNEL 

005F 

1 

016C 

DC 

AREAG 

0060 

0 

0000 

r.'PON 

DC 

*  -  * 

0001 

0 

CO  F  2 

LD 

ADO*  1 

RESTORES  ADR.  II.  AREAG 

0002 

01 

94000238 

S 

L 

pouts 

0004 

01 

1)4000054 

STO 

L 

ADO*  1 

0060 

0 

C0E4 

LD 

LC0P1*1 

0007 

01 

E400022E 

A 

L 

TfJPTS 

0009 

0 

00E1 

STO 

L00P1 ♦ 1 

OOOA 

20 

26063A15 

SKI  PP 

ligf 

VJRTYI! 

CONVERSION  CHANNEL  3 

0000 

0 

2001 

DC 

/  2001 

006C 

1 

0488 

DC 

CARIG 

0  0  00 

0 

0000 

DC 

0 

0006 

20 

2G6C3A15 

LI  CF 

WRTYN 

0  0  OF 

0 

0001 

DC 

7  0001 

0070 

0 

7  0  F  0  . 

IIDX 

*-  3 

00  71 

20 

2GG63A15 

LI  OF 

WRTYN 

HEAOING  FOR  CHANNEL  3 

0072 

0 

2001 

DC 

/  2001 

0073 

1 

04A2 

DC 

DMES2-1 

0074 

0 

0000 

DC 

0 

00  75 

20 

2GG63A15 

L  1  OF 

WRTYII 

0070 

0 

0001 

DC 

/0001 

0077 

0 

70  FD 

MOX 

»-3 

0078 

20 

2GGG3A15 

I.IBF 

WRTYN 

CARRAGE  RETURN 

0070 

0 

2001 

DC 

/  2 001 

00  7 A 

1 

0488 

DC 

CARIG 

007P 

0 

0000 

DC 

0 

007C 

20 

26603A15 

LIBF 

WRTYN 

0070 

0 

0001 

DC 

/0001 

007E 

0 

70  FD 

MDX 

•  -3 

007F 

01 

G5S0022E. 

I.DX 

1  1 

TNPTS 

0081 

01 

C4000485 

LD 

L 

AORO 

SETS  ADR.  FOR  3. CHANNEL 

0083  01  0400022F 

S 

L 

FNPTS 

00C5 

01 

94000486 

S 

L 

DEC1 

0087 

0 

0001 

STO 

LOOP2+1 

0088 

00 

CCOOOOOO 

L00P2 

LDD 

L 

*  -  * 

LOOP  FOR  CHANNEL  3 

008A 

30 

0  5 1 0  G  4  E  3 

CALL 

EDFLT 

FLOAT  THE  VALUE 

008C 

20 

05109940 

LIBF 

EOIV 

DIVISION  ROUTINE 

QOSP 

1 

0230 

DC 

DVSER. 

008E 

30 

05898G40 

CALL 

ESQR 

SQUARE  ROOT 

0090 

20 

091S99C0 

LIBF 

1  F  IX 

CONVERTS  TO  INTEGER 

0091 

01 

nuoooicc 

ADI 

STO 

L 

AREAII 

FIELD  FO  1  NT .  VALUE 

0003 

01 

74010092 

MDX 

L 

AD1*1,1 

0095 

01 

74020089 

MDX 

L 

L00P2+1, 2 

0097 

0 

7102 

100PB 

MDX 

1 

2 

113 


0098 

0 

70EF 

MDX 

L00P2 

f)Qqq 

01 

ruQGn23fi 

in 

L 

PONTS  . 

oo  n". 

6 

0002 

STO 

tIPPOtl 

oooc 

20 

23A17171 

LI3F 

TYPE1 

PRINTOUT  3. CHANNEL 

009D 

1 

01CC 

nc 

AREAH 

C09E 

0 

0000 

IIPPON 

00 

*  -  * 

OOOF 

0 

COF  2 

AD1  +  1 

RESTORES  l.AOR.  IN  ARfeAH 

OOAO 

01 

34000238 

s 

L 

pcnts 

OOA2 

0 

DOEF 

ADI  ♦  1 

00A3 

01 

OCOOOl 1 2 

Itl 

L 

IOCCA 

READ  DATA  ENTRY  SWITCH 

0  0/.  3 

01 

E4000114 

L 

MAS  C 

OOA  7 

4808 

4**i 

■f 

0  OA  8 

0 

7006 

MDX 

OUT 

00A9 

20 

05059115 

LI  OF 

CARD!) 

FEEDS  1  CARD 

OOAA 

0 

3000 

DC 

/  3000 

00A3 

0 

0000 

DC 

0 

OOAC 

20 

03059115 

L  l  b  F 

CAROM  4 

0 

COOP 

DC 

/oooo 

OOAE 

0 

70FU 

MDX 

•-3 

OOAF 

01 

C4Q0023A 

ID 

L 

CLICK 

PUNCH  RUNNO 

0031 

20 

02255103 

LI  BF 

G  IUDC 

0  00  2 

1 

0117 

DC 

OTPUT 

00"3 

01 

C4000233 

LO 

L 

WAVES 

003  3 

20 

LI  OF 

G  IIJDC 

OOOC 

1 

011D 

DC 

OTPUT+6 

000  7 

01 

C4000169 

LD 

L 

DEC12 

ooro 

0 

00  5  C 

STO 

AREA 

OOCA 

20 

03059115 

LIDF 

CAR  DM 

OOBB 

0 

2000 

DC 

/2000 

1 

0116 

DC 

AREA 

0  000 

0 

0000 

DC 

0 

QOT  F. 

20 

03059115 

LI  OF 

CARDJI 

00nF 

0 

0000 

DC 

/OOOO 

QQCO 

0 

70FD. 

MDX 

*-  3 

00C1 

01 

04 0001 6A 

LD 

L 

DEC72 

PUNCH  MEAN  AND 

0 

l»7*WWI 

STO 

AREA 

PERIODIC  ON  CARDS 

oo:u 

01 

C5S00238 

LDX 

11 

POHTS 

0006 

0 

620C 

PUf.'Ctl 

LDX 

2 

12 

12  VALUES  PER  CARD 

0007 

01 

C40001CC 

COMV 

LD 

L 

AREAII 

00C3 

20 

02255103 

LI  OF 

DINDC 

OOCA 

1 

0117 

OTP 

DC 

OTPUT 

01 

MDX 

L 

OTP,  6 

FIELD  FOR  NEXT  VALUE 

OOCD 

01 

740100C8 

MDX 

L 

couv+i, l 

NEXT  VALUE 

OOCF 

0 

72FF 

MDX 

2 

-l 

OODO 

0 

70F6 

MDX 

COIIV 

00D1 

20 

03059115 

IIBF 

CARDN 

00D2 

0 

2000 

DC 

/  2000 

■mi 

1 

0116 

AREA 

0004 

0 

0000 

DC 

0 

00D5 

20 

03059115 

LI  BF 

CARDN 

0006 

0 

0000 

DC 

/OOOO 

0007 

0 

70FD 

MDX 

.-3 

0008 

01 

74C300CA 

MDX 

L 

OTP, -72 

RESTORES  1.  OUTPUT  AD. 

OODA 

0 

71F4 

MDX 

1 

-12 

OODP 

0 

70  FA 

MDX 

PUNCH 

.  OODC 

:  oi 

C4000238 

LO 

L 

POUTS 

PRINTOUT  OF  REST 

QODE 

0 

1890 

SRT 

16 

OOOF 

01 

ACC00169 

D 

L 

DEC12 

00E1 

0 

1090 

SLT 

16 

00F2 

oi_ 

B  4  0  0  0 1 6  7 

C".P 

I.  . 

ZERO 

00F4 

0 

7002 

MDX 

GOON 

00F5 

o 

7023 

MDX 

P.STOR 

00E6 

0 

7022 

MDX 

RSTOR 

OOF.  7 

01 

D400016B 

GOON 

STO 

L 

REST 

00E9 

01 

A  4 000168 

M 

L 

0EC6 

OOFB 

0 

1090 

SLT 

16 

114 


OOEC 

01 

D 4 000116 

STO 

1 

AREA 

88P 

i  n 

SET  ADRESS  IN  AREAH 

8 

8888 

STO 

CON  VV* 1 

00  F  0 

01 

64000160 

A 

L 

REST 

00  F  2 

0 

D005 

STO 

CON  V* 1 

OOF  0 

01 

66800160 

LDX 

12 

REST 

00F5 

00 

C4000000 

convv 

LD 

L 

*  -  # 

OOF7 

20 

02255103 

L 1  BF 

6  IUDC 

OOFS 

1 

0117 

OTTP 

DC 

OTPUT 

OOFO 

01 

74Q600F8 

MDX 

L 

OTTP, 6 

00  FO 

01 

740100F6 

MDX 

L 

C0NVV+1, 1 

OOFP 

0 

72FF 

MDX 

2 

-1 

OOFu 

0 

70  F  6 

MDX 

CONVV 

OOFF 

20 

03059115 

LI6F 

CAR  Dll 

0100 

0 

2000 

DC 

/  2  000 

0101 

1 

0116 

DC 

AREA  _ 

0102 

0 

0000 

DC 

0 

oim 

20 

03050115 

LI  CF 

CARDN 

nm 

0 

0000 

DC 

/OOOO 

0105 

0 

70  FD 

MDX 

*-3 

0106 

0 

C0F1 

LD 

OTTP 

RESTORES  OTTP 

0107 

0 

900F 

S 

AREA 

OIOS 

0 

DOFF 

STO 

OTTP 

0109 

01 

C40000C8 

R5T0R 

LD 

L 

CCNV+1 

RESTORES  ADR.  IN  CONV+l 

0106 

01 

94000238 

r 

O 

L 

POMTS 

010D 

01 

D40000CS 

STO 

L 

CONV+1 

010F 

30 

059C‘J8C0 

OUT 

CALL 

EXIT 

0112 

0000 

nss 

E 

0 

0112 

0 

0000 

I  OCCA 

DC 

/  cooo 

READ  DATA  ENTRY  SW . 

0113 

0 

0740 

DC 

/  0  7  4  0 

0114 

0 

0001 

MAS  C 

DC 

/  0001 

0115 

0 

0002 

;;asc2 

DC 

/  000  2 

MASK  FCR  DES  2 

0116 

c 

0000 

AREA 

DC 

*  -  * 

output  area  for  cards 

0117 

0050 

OT  PUT 

BSS 

80 

0167 

0 

0000 

ZERO 

DC 

0 

01CE 

0 

0006 

DECO 

DC 

6 

0169 

0 

OOOC 

DEC12 

DC 

12 

016A 

0 

0048 

DEC72 

DC 

72 

016B 

0 

0000 

REST 

OC 

*  -  * 

016C 

0060 

ARCAC 

CSS 

96 

STOR.  1.  CHANNEL 

010C 

0060 

AREAII 

CSS 

96 

STOR.  2.  CHANNEL 

0  2  2C 

0 

0009 

DEC9 

DC 

9 

0  2  2D 

0 

0000 

MNPTS 

DC 

*  -  * 

0  22F 

0 

0000 

TNPTS 

DC 

*-  * 

022F 

0 

0000 

FfIPTS 

DC 

*  -  * 

02  30 

0003 

DVSCR 

RSS 

3 

0234 

0000 

BSS 

E 

0 

0254 

31 

C6253COO, 

1  NFO.” 

dsa 

F  1  L20 

INFORMAT IONF  FOR  EXP. CONST 

0237 

0 

0000 

DELAY, 

,oc 

*  -  * 

0  2  3  S 

0 

0000 

POUTS 

DC 

*  -  # 

0  239 

0 

0000 

l.'AVFS 

DC 

0  2  3A 

0 

0000 

CLICK 

DC 

*  -  « 

0236 

0006 

CSS 

6 

0242 

0000 

css 

_E _ 

0 

0242 

31 

06253CB  2 

CPVST 

DSA 

F  1  L  2  2 

USES  FIL22  AND  FIL23 

0  245 

0240 

AREAD 

css 

576 

SQUARE  OF  RANDOM  FLUCTU. 

0485 

1 

0245 

ADRU 

DC 

AREAD 

0486 

0 

0001 

DEC1 

DC 

1 

0487 

0 

0240 

wor.oc 

DC 

576 

_qiaS£  n 

nn  m _ 

CAR  10 

DC 

1 

0  489 

0 

8181 

DC 

/  81 81 

043A 

0 

0008 

DC 

8 

0480 

0008 

OUTPT 

CSS 

8 

0493 

0 

2121 

DC 

/  2121 

0494 

0 

OOOD 

DC 

OMESA-PMES 

1 

0495 

0018 

DME31 

DMF.S 

’ 2RRMS  VALUES  OF  PRESSURE' E 

■fA 


»  '  «|y 


